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INTRODUCTION 

A  prevalent  practice  in  beef  packing  plants  is  to  slaughter 
the  animal   and  chill   the  carcass  for  approximately  24-72  hours.     Fre- 
quently, following  the  cooling  period,  the  excess  fat  and  bone  are 
removed  as  the  carcass  is  fabricated  into  muscles  and  muscle  systems 
before  packaging  and  shipment.     Recent  studies  on  product  quality 
indicate  that  a  procedure  usually  termed  hot  boning  or  hot  processing 
can  be  applied  to  produce  beef  steak,   roast,  and  ground  products. 
Hot  processing  entails  the  removal   of  excess  fat  and  bone  from  the 
beef  carcass  soon  after  slaughter.     The  traditional   cooling  period 
after  slaughter  is  either  shortened  or  eliminated,  thus,  the  carcass 
is  warm  when  cut  and  the  resultant  cuts  are  subsequently  chilled. 

This  study  compares  two  hot-processing  procedures  to  conven- 
tional  centralized  beef  processing  to  determine  what  method  is  the 
most  energy  efficient.     This  comparison  determines  the  difference  in 
the  average  energy  requirements  of  the  various  options  and  is  not 
applicable  to  equipment  sizing. 

The  energy  comparisons  are  based  upon  model   facilities  that 
slaughter  sixty  head  of  cattle  an  hour  for  eight  hours  a  day.     How- 
ever, conclusions  found  for  this  size  hot-processing  facility  should 
not  be  unreservedly  ascribed  to  other  sized  plants. 
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OBJECTIVES 

1)  To  determine  the  rate  of  product  cooling  and  the  internal  product 
temperature  of  the  proposed  hot-processing  schemes  through  transient 
heat  analysis. 

2)  To  develop  a  computer  program  that  calculates  the  energy  requirements 
of  conventional  and  hot-processing  techniques. 

3)  To  compare  the  energy  requirements  of  two  different  hoi-processing 
schemes  against  conventional   cold  processing. 


REVIEW  OF  LITERATURE 

Conventional   and  Hot  Processing 

TyDically,   in  a  packing  plant  fabricating  boxed  beef  products, 
immediately  after  slaughter  the  carcass  is  placed  in  a  chill  or  drip 
cooler  maintained  near  30°F  for  18-24  hours.     At  the  expiration  of 
this  period  the  carcass  is  moved  to  a  holding  cooler  operating  near 
32°F  where  it  may  remain  up  to  72  hours  after  slaughter.     The  carcass 
is  then  moved  to  the  boning  line  where  excess  fat  and  bone  are  removed 
and  the  carcass  is  fabricated  into  the  desired  cuts.     The  boning  line 
temperature  is  approximately  50°F  and  the  time  the  carcass  spends  on 
the  line  is  somewhat  variable.     The  fabricated  product  is  then  vacuum- 
packed.     It  is  placed  in  boxes  and  stored  in  a  28-32°F  cooler  awaiting 
shipment  to  retail  markets.     Throughout  this  study  this  conventional 
method  of  processing  will   be  termed  cold  processing. 

Henrickson   (1975)  stated  that  the  refrigerated  space  required 
for  chilling  a  600  pound  beef  carcass  is  86,400  cu  in.     He  added  34,000 
cu  in    of  space  above  and  below  the  carcass  making  a  total   of  120,000 
cu  in.     He  assarted  that  the  edible  portion  of  the  carcass  could  be 
cooled  in  26,000  cu  in  of  space  resulting  in  a  78  percent  savings  in 
space.     A  30  percent  reduction  in   the  energy  requirement  to  cool   each 
carcass  was  also  claimed  by  Henrickson  on  the  basis  that  22.380  Btu 
will   reduce  the  edible  portion  of  the  carcass  from  102°F  to  32°F,  while 
31,500  Btu  would  be  required  to  achieve  the  same  temperature  reduction 
in  the    total   carcass. 

These  estimated  energy  savings  coupled  with  a  recent  and  growing 
conservation  necessity  had  led  to  intensified  interest  in  hot  processing. 
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The  intent  of  hot  processing  is  to  remove  excess  fat  and  bone  as  soon  as 
possible  following  slaughter  without  adversely  affecting  product  quality. 
This  would  eliminate  the  use  of  energy  to  chill   carcass  portions  which  will 
be  rendered. 

A  variety  of  hot-processing  schemes  have  been  proposed  and  the  primary 
thrust  of  the  research  deals  with  the  quality  aspect  of  these  proposals.     Two 
basic  approaches  appear  feasible.     One  approacn  is  to  condition  the  carcass 
at  6C°F  for  selected  periods  of  up  to  10  hours  post-mortem  and  to  then  hot 
process  the  sides.     Research  has  shown  that  sides  held  intact  from  3-8  hours 
post-mortem,  then  fabricated,  give  cuts  equal   to  the  conventional  treatment 
considering  yield,  tenderness,  color  and  flavor  (Kastner  et  al ,   1973;  Kastner 
and  Russell,  1975;  Kastner,   1976;  Will  and  Henrickson,   1976). 

The  second  basic  approach  is  to  cut  the  carcass  soon  after  slaughter 
and  to  subsequently  condition  the  carcass  components.     This  technique  nor- 
mally results  in  acceptable  organoleptic,  yield  and  bacterial   characteristics 
(Schmidt  and  Gilbert,  1970  and  Schmidt  and  Keman,  1974).     A  variation  of  this 
technique  which  is  being  widely  considered  is  electrical   stimulation  of  the 
carcass  to  hasten  rigor  mortis  and,  thus,  allow  faster  chilling  and  excision 
(McCollum  and  Henrickson,  1977;  Gilbert  and  Davy,  1976). 

One  of  the  still   largely  unanswered  questions  of  hot  processing  is 
the  potential  bacterial  problem.     Because  hot  processing  involves  vacuum 
packaging  hot  cuts  and  immediately  boxing  and  chilling  them  the  effect  of 
this  increased  product  temperature  on  the  microbial   count  must  be  investigated 
Schmidt  and  Gilbert  (1970)  concluded  that  pre-rigor  excision   (1-2  hours  post- 
mortem) with  subsequent  aging  in  vacuum  packages  at  59°F  for  24  hours  could 


produce  both  organoleptically  and  microbiologically  acceptable  beef.  Fung 
et  al  ( 1979)  found  that  meat  that  was  hot  processed  and  stored  for  14  days 
was  considered  bacteriologically  acceptable. 

Cooling  Load  Analysis 

Determination  of  the  total  amount  of  energy  required  to  chill  the 
product  and  air  condition  the  necessary  physical  plant  requires  considera- 
tion of  seven  major  areas  of  heat  gain.  Two  are  inherent  to  the  structures, 
heat  conduction  through  exterior  walls  and  roof,  and  heat  conduction  through 
interior  partitions.  Ventilation  and  infiltration  constitute  another  area 
of  heat  gain.  Heat  will  also  be  acquired  within  the  structure  due  to  people, 
lights  and  other  electrical  equipment.  Finally,  the  refrigeration  energy 
required  to  cool  the  product  must  be  added  to  these  sources  of  heat  gain. 

Traditional  techniques  (ASHRAE  Fundamentals,  1977)  of  evaluating 
cooling  loads  were  generated  for  the  purpose  of  equipment  sizing.  These 
were  based  on  so-called  "design  conditions."  This  is  the  worst  possible 
situation  which  could  ever  be  placed  on  the  system.  Since  these  condi- 
tions would  only  exist  for  a  short  time,  if  at  all,  it  exaggerates  the 
amount  of  refrigeration  energy  actually  required  by  the  system. 

Through  the  use  of  digital  computers  it  is  possible  to  rigorously 
calculate  energy  requirements  by  solving  heat  balance  equations  at  all 
interior  surfaces  of  a  given  space  (Mitalas  and  Stephenson,  1967).  This 
procedure  provides  the  best  method  of  determining  the  energy  requirements 
of  each  unique  system  but  involves  a  time  consuming  and  laborious  technique. 


Kusuda  (1976)  presents  one  of  the  few  computer  programs  that  rigorously 
determines  instantaneous  cooling  load. 

The  transfer  function  concept  of  determining  cooling  requirements 
was  introduced  by  Mitalas  and  Stephenson  (1967)  as  a  simplification  of 
their  rigorous  method  described  above.  This  procedure  was  generated  by 
calculating  the  cooling  load  by  the  rigorous  method  for  typical  construc- 
tions. The  components  such  as  solar  heat,  conduction  heat  gain,  lighting, 
equipment  and  personnel  load  are  simulated  by  pulses  of  unit  strength. 
Transfer  functions  are  then  calculated  as  numerical  constants  which  represent 
the  cooling  load  corresponding  to  the  input  exaltation  pulse.  Transfer 
functions  are  generated  for  a  number  of  typical  constructions  and  then  are 
assumed  to  be  independent  of  the  input  pulses.  The  transfer  function  can  now 
be  multiplied  by  a  time  series  representation  of  heat  gain  to  obtain  the 
cooling  load  (ASHRAE  Fundamentals,  1977).  The  transfer  function  technique 
provides  an  answer  of  reasonable  accuracy  and  embodies  a  simpler  approach 
than  the  rigorous  method. 

ASHRAE  Fundamentals  (1977)  specifies  the  accepted  standard  for  cal- 
culating the  heat  gain  due  to  personnel,  product,  lighting,  equipment, 
ventilation  and  infiltration.  These  methods  will  be  used  and  discussed 
in  the  thermal  analysis  section. 


THERMAL  ANALYSIS 

Hot-Processing  Option  Overview 

Two  hot-processing  approaches  were  selected  for  analysis.     These 
schemes  reflect  the  two  basic  hot-processing  procedures  described  previously. 
Figure  1   shows  the  first  option  where  the  carcass  is  condi toned  for  8  hours 
and  then  fabricated.     The  second  hot  processing  option  is  shown  in  Figure 
2.     This  option  involves  cutting  the  carcass  soon  after  slaughter  with  a 
subsequent  conditioning  and  chill  period.     Appencices  A  and  3   show  the 
calculation  of  the  cooling  rates  of  the  product  for  each  option.     Under 
hot-processing  option  I  the  internal   temperature  would  reach  28.3CF  in 
104  hours  and  under  option   II   it  would  be  28.3  °F  in  96  hours. 

Figure  3  depicts  the  conventional   cold-processing  method.     An 
analysis  of  its  cooling  rate  is  evaluated  in  Appendix  C.     It  requires 
approximately  72  hours  for  the  internal  product  temperature  to  reach 
28.3°F  under  this  procedure. 

A  comparison  of  the  energy  requirements  of  these  three  different 
schemes  of  meat  packing  entails  evaluating  the  components  that  constitute 
a  source  of  heat  gain  for  the  system.     This  analysis  does  not  attempt  to 
evaluate  the  total  amount  of  energy  used  for  each  packing  method.     The 
total  quantity  of  energy  used  is  dependent  upon  each  unique  building  design 
and  does  not  lend  itself  to  the  generalized  approach  attempted  here.     Rather, 
the  objective  of  this  research  is  to  compare  areas  of  dissimilarity  within 
each  system  and  evaluate  the  difference  in  the  energy  requirements  of  each 
system.     The  slaughtering  operation,  boning  line,   and  packaging  and  boxing 
phase  will  not  be  considered  as  these  are  assumed  identical   for  each  oper- 
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tion.  Carcass  fabrication  and  meat  packaging  occur  so  rapidly  that  heat 
transfer  during  these  operations  is  virtually  nonexistent. 

Appendix  D  contains  a  sample  calculation  of  each  component  of  the 
refrigeration  energy  load. 
Product  Load 

The  refrigeration  load  imposed  by  the  product  results  from  the 
energy  necessitated  to  reduce  the  product  temperature  from  receiving  to 
leaving  temperature.  This  load  can  be  evaluated  usinq  equation  1, 

q  =  mc(trt2)  (!) 

where  q  is  the  quantity  of  heat  removed,  m  is  the  mass  of  the  product, 
c  is  the  specific  heat  of  the  product  above  freezing,  and  t*   and  t«, 
respectively,  are  the  temperatures  cf  the  product  as  it  enters  and  leaves 
the  cooler.     This  equation  is  valid  as  long  as  the  meat  doesn't"  reach  the 
freezing  temperature  of  28°F  (ASHRAE  Fundamentals,  1977d). 

ASHRAE  Applications   (1974a)  gives  typical   temperature  cooling  curves 
for  the  carcass  in  the  chilling  and  holding  coolers  of  the  cold-processing 
method.     However,  no  similar  curves  are  available  for  the  carcass  as  it  is 
held  in  tne  higher  temperature  conditioning  room  characteristic  of  hot- 
processing  option  I.     Therefore,  a  transient  heat  analysis  of  the  carcass 
under  all   conditions   is  necessary  to  arrive  at  the  average  product  tempera- 
ture under  the  same  methods.     Appendices  A,  3  and  C  show  this  calculation 
for  hot-processing  option  I,  hot-processing  option  II  and  cold-processing, 
respectively. 

The  mass  of  the  carcass  will   decrease  once  it  reaches  the  oonino 


12 
line  due  to  removal  of  the  excess  fat  and  bone.  A  typical  choice  beef 
carcass  consists,  on  a  weight  basis,  of  51  percent  lean,  15  percent  bone 
and  34  percent  fat  (Watt  and  Merrill,  1963).  For  this  study  it  was  assumed 
that  excess  fat  (12  percent)  and  all  bone  would  be  eliminated  once  the 
carcass  completed  the  boning  line  process.  Under  all  schemes  (hot  and 
cold)  it  is  proposed  that  the  deboned  carcass  be  cut  into  the  desired  sub- 
primals  and  packaged  in  vacuum  bags.  The  packaged  cuts  are  then  placed  in 
cardboard  boxes  and  stored  in  the  shipping  cooler. 

The  specific  heat  of  the  meat  varies  due  to  the  change  in  the  carcass 
percentages  of  lean,  fat  and  bone.  It  was  calculated  that  after  boning,  the 
boxed  cuts  from  the  carcass  averaged  70  percent  lean  and  30  percent  fat. 
Tne  new  specific  heat  using  a  value  of  0.50  Btu  per  (lb)  (F  deg)  for  fat 
and  C.80  Btu  per  (lb)  (F  deg)  for  lean  is  0.71  Btu  per  (lb)  (F  deg)  com- 
pared to  the  original  carcass  specific  heat  of  0.75  Btu  per  (lb)  (F  deg) 
(ASHRAE  Applications,  1974a). 


Shrouding 

A  carcass  that  is  conventionally  cold  processed  is  usually  shrouded 
to  enhance  appearance  by  smoothing  and  bleaching  the  surface  fat.     A 
shroud  is  a  cloth  saturated  in  120-130CF  weak  brine.     It  remains  on  the 
carcass  throughout  the  initial  chilling  period  but  is  removed  before  the 
carcass  is  transferred  to  the  holding  room.     Shrouds  lose  about  six  pounds 
of  water  per  carcass  to  the  chill   cooler,  largely  due  to  evaporation.     To 
determine  the  heat  gain  in  the  chill   room  it  is  necessary  to  consider  both 
the  sensible  and  latent  loads  imposed  by  the  warm,  wet  shroud. 

The  sensible  load  can  be  determined  by  equation  2, 

Sensible  =  mwcw  (ts  "  V  (2) 

where  q  is   the  amount  of  heat  removed,  m,    is   the  amount  of  water  loss,  c 

w  w 
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is  the  specific  heat  of  water,  t  is  the  temperature  of  the  shroud  water 
and  t  is  the' temperature  of  the  chill  cooler. 

The  latent  load  imparted  by  the  shroud  is  calculated  by  equation  3, 

qlatent  =  hfg  {mJ  (3) 

where  q  is  the  amount  of  heat  to  be  removed,  h^  is  the  heat  of  vaporiza- 
tion of  the  water  and  m  is  the  amount  of  water  lost  by  the  shroud.  An 

w 

average  value  of  h^  between  the  125°F  temperature  the  shroud  enters  the 
chill  cooler  and  the  30°F  temperature  the  shroud  leaves  the  cooler  is 
1049.5  Btu  per  (lb)  water. 

Building  Transmission 

The  conduction  losses  through  the  roof  and  exterior  walls  will  be 
determined  by  the  previously  described  transfer  function  concept.  The 
heat  gain  to  the  space  from  the  exterior  walls  and  roof  is  determined  by 
equation  4,  (ASHRAE  Fundamentals,  1977), 

"e.r  ■  A  [n:o  >n  <W>  '  £  \  !W  -  Vc  £  «     (4) 

where  q         is  the  hourly  heat  gain  to  the  conditioned  space;  A  is  the 

exposed  area  of  the  roof  or  wall;  r  is  the  time  the  calculation  is  being 

done;  A  is  the  time  interval   between  calculations;  n  is  the  summation 

index;  t      is  the  constant  indoor  air  temperature;  b   ,  c     and  d     are  the 
rt-  n       n  n 

transfer  function  coefficients;  and  t  A   is  the  sol-air  temperature 

e,r-nA  r 

at  time  r-nA. 

The  building  transmission  calculation  will  be  done  the  twenty-first 
day  of  each  month  at  one  hour  intervals  beginning  at  1:00  A.M.     Thus  the 
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r  value  will  proceed  from  1:00  A.M.  to  2:00  A.M.,  etc.  with  a  A  value  of 
one  hour.  The  number  of  summation  indices  (n)  in  this  calculation  is  a 
function  of  the  building  construction.  ASHRAE  Fundamentals  (1977)  gives 
this  information  based  upon  many  common  industrial  constructions.  The 
values  of  the  b  ,  d  and  c  terms  can  be  read  from  the  same  source  once 
the  desired  type  of  building  construction  is  selected. 

The  (q  r_nA)  term  is  the  value  of  the  heat  gained  from  the  calcu- 

lation  performed  for  the  previous  hour.  Initially  this  value  is  set  at 
zero  and  the  calculation  is  then  iterated  sufficiently  to  negate  the  effect 
of  this  procedure. 

The  value  of  (t     )   or  the  sol-air  temperature  one  hour  previous 
to  when  the  calculation  is  performed  is  calculated  by  equation  5,  (ASHRAE 
Applications,  1974b), 

aIt    EAR  ,,v 

v^-—  "V  <5) 

where  the  sol -air  temperature  is  t  ,  t  is  the  outdoor  air  temperature,  a  is 
the  absorptance  of  the  surface  for  solar  radiation,  I  is  the  total  solar 
radiation  incident  upon  the  surface,  h  is  the  coefficient  of  heat  transfer 

by  long  wave  radiation  and  convection  at  the  outer  service,  E  is  the 

hemispherical  emittance  of  the  surface  and  AR  is  the  difference  between 

the  longwave  radiation  incident  on  the  surface  from  the  sky  and  surroundings, 

and  the  radiation  emitted  by  a  blackbody  at  outdoor  air  temperature. 

The  value  of  t  should  be  gathered  from  historical  weather  data  for 
o 

an  area  near  the  plant  location.  These  should  be  the  1:00  A.M.,  2:00  A.M., 
etc.  temperature  readings  to  correspond  to  the  previously  selected  values 

of  r.  The  values  of  -x-  given  in  ASHRAE  Fundamentals  (1977)  are  0.15  for 

o 
a  light-colored  surface  and  0.30  for  a  dark-colored  surface.  Any  interme- 
diate colored  surface  should  be  interpolated  between  these  values.  I.  is 
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given  in  ASHRAE  Fundamentals  (1977)  as  equal  to  1.15  multiplied  by  the 
solar  heat  gain  factor  (SHGF).  The  values  of  the  SHGF  are  given  in 
ASHRAE  Fundamentals  (1977b)  as  a  function  of  latitude,  date,  solar  time 
and  building  wall  orientation.  Solar  time  may  run  either  faster  or 
slower  than  civil  time  and  this  difference  can  be  significant.  Its  cal- 
culation is  described  in  Appendix  E.  Table  1  in  Appendix  E  lists  the 
"equation  of  time"  variable  necessary  to  calculate  solar  time  while  Table 
2  shows  the  solar  time  throughout  the  year  for  Manhattan,  Kansas. 

The  hemispherical  emittance,  E,  is  given  in  ASHRAE  Fundamentals 
(1977a)  based  upon  the  building  surface.  Values  of  AR  for  horizontal 
surfaces  that  receive  longwave  radiation  from  the  sky  only  are  given  as 
20  Btu  per  (hr)  (sq  ft)  and  for  a  vertical  surface  AR  is  considered  to 
be  zero  (ASHRAE  Fundamentals,  1977). 

The  heat  load  which  will  be  imposed  due  to  transmission  from  in- 
terior partitions,  ceiling  and  floors  is  given  by  equation  6, 

«D.r  =  UA  (tb  "  Vc>  (6) 

where  q    is  the  hourly  heat  gain  through  the  interior  partitions,  U  is  the 
p»r 

coefficient  of  overall  heat  transfer  between  the  adjacent  space  and  the 
cooler  space,  A  is  the  exposed  area,  t,  is  the  temperature  of  the  adjacent 
space  and  t   is  the  temperature  of  the  cooler  space. 

The  U  value  is  dependent  upon  the  construction  type.  It  may  be  read 
directly  in  ASHRAE  Fundamentals  (1977)  based  upon  a  construction  descrip- 
tion or  it  may  be  determined  by  summing  the  U  values  of  the  individual 
construction  components  as  described  in  ASHRAE  Fundamentals  (1977a). 

The  heat  transfer  from  floors  in  direct  contact  with  the  ground, 
or  over  an  underground  basement  that  is  neither  ventilated  nor  warmed, 
may  be  neglected  for  cooling  load  estimates  (ASHRAE  Fundamentals,  1977). 
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Infiltration 

Infiltration  is  air  leakage  through  building  cracks,  around  windows 
and  doors,  and  through  building  floors  and  walls.  The  amount  of  infiltra- 
tion depends  upon  the  type  of  building  construction,  workmanship,  and 
condition  of  the  building.  Air  leakage,  in  refrigerated  rooms,  is  pri- 
marily due  to  opening  doors.  Each  time  a  door  is  opened  outside  air  enters 
the  cooled  space.  The  temperature  and  moisture  content  of  this  infiltrating 
air  must  be  brought  to  the  refrigerated  space  condition,  which  constitutes 
a  sensible  and  latent  load  to  the  system.  The  volume  of  air  entering  a 
walk-in  refrigerator  due  to  door  openings  is  estimated  by  ASHRAE  Fundamentals 
(1977c)  based  upon  practical  field  experience. 

The  total  load  component  resulting  from  air  infiltration  may  be 
calculated  by  equation  7, 

- 

q  =  4.5   (cfm)    (Ah)  (7) 

where  q  is  the  hourly  heat  load,  Ah  is  the  difference  in  enthalphy  between 
inside  and  outside  air,  and  cfm  is   the  volume  of  air  entering  the  space,   ft/min. 
The  4.5  factor  is  strictly  valid  only  for  air  with  a  density  of  0.075  lb 
dry  air  per  cubic  foot.     However,  accuracy  is  not  noticably  affected  by 
its  use  throughout. 

Enthalphy  values  for  the  air  within  the  refrigerated  space  are  shown 
in  Table  3  for  the  coolers  considered  in  this  study.     The  enthalphy  of  the 
adjoining  office  space  is  23.7  Btu  per  pound  of  dry  air,  assuming  room  con- 
ditions of  70°F  and  40%  relative  humidity.     Enthalphy  values  are  a  function 
of  temperature  and  relative  humidity.     Outside  enthalphy  values  are  cal- 
culated by  the  computer  code  described  later  by  inputting  the  dry  bulb  air 
temperature  and  relative  humidity  values  at  the  time  desired. 
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Table  3.  Enthalphy  Values  Used  in  Infiltration 
Calculations  -  Coolers 


Cooler  Temp.  Cooler  Relative  Humidity      Cooler  Enthalphy 
(F)               (%)  (Btu  per  lb  dry  air) 

28  85  9.6 

30  85  10.3 

32  85  11.1 

60  80  24.0 
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Personnel  Load 


Humans  give  off  moisture  and  heat  that  subsequently  constitutes 
a  load  to  the  system.   The  rate  that  heat  and  moisture  is  given  off 
depends  upon  space  temperature,  type  of  work  or  activity,  clothing,  size, 
etc.  ASHRAE  Fundamentals  (1977c)  gives  the  average  hourly  system  load 
due  to  occupancy  as  a  function  of  cooler  temperature.  The  heat  equivalency, 
as  given  by  this  source,  for  a  person  in  a  60  F  cooler  is  600  Btu  per  hour 
(by  extrapolation)  and  within  a  30°F  cooler  the  heat  equivalency  per  person 
is  950  Btu  per  hour.  These  values  should  be  increased  for  a  heavy  traffic 
load  where  many  people  may  enter  the  refrigerated  space  for  a  short  time 
and  carry  with  them  heat  above  that  indicated. 

Electrical  Load 

The  two  major  contributors  to  electrical  load  considered  here  are 
lights  and  motors.  The  electrical  energy  dissipated  within  the  cooled 
space  due  to  motors  will  depend  upon  the  size  (horsepower  rating)  of  the 
motor  and  the  placement  of  the  motor.  ASHRAE  Fundamentals  (1977c)  gives 
the  heat  equivalent  of  electrical  motors  where  both  the  useful  output  and 
motor  losses  are  dissipated  within  the  refrigerated  space;  or  where  the 
motor  losses  are  dissipated  outside  refrigerated  space  and  useful  work  is 
expended  within  refrigerated  space;  or  where  motor  losses  are  dissipated 
within  the  refrigerated  space  and  useful  work  is  expended  outside  of  the 
refrigerated  space.  It  will  be  assumed,  for  this  case,  that  motor  losses 
are  within  the  refrigerated  space  and  useful  work  is  expended  within 
the  refrigerated  space.  The  assumed  fan  horsepower  for  this  example  will 
be  36  for  the  chill  and  shipping  cooler.  The  holding  cooler  and  conditioning 
room  will  be  assumed  to  have  a  24  horsepower  rating.  The  hourly  heat  load 
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due  to  the  electrical  motor  is  the  product  of  the  fan  horsepower  rating 
and  the  conversion  factor  to  the  desired  energy  units,  in  this  case 
2545  Btu  per  (hp)  (hr). 

The  energy  load  to  the  refrigerated  space  as  a  result  of  electri- 
cal lighting  can  be  determined  by  equation  8. 

q  =  total  wattage  x  use  factor  x  special  allowance  factor  x  3.41    (8) 

The  use  factor  is  the  ratio  of  the  wattage  in  use  to  the  total  installed 
wattage.   The  special  allowance  factor  is  used  for  fixtures  that  actually 
require  more  energy  than  their  rated  wattage,  such  as  fluorescent  or  rapid 
start  fixtures.  The  3.41  figure  is  the  conversion  factor  from  watts  to 
Btu  per  hour. 
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DESCRIPTION  OF  PROGRAM 

A  computer  program  called  HEATLD  was  written  to  calculate  the  energy 
requirements  of  each  processing  option.  It  contains  a  main  program,  three 
subroutines  which  determine  the  energy  consumption  of  each  of  the  three 
processing  options,  and  a  fourth  subroutine  (TRANEX)  which  calculates 
the  external  transmission  using  the  transfer  function  coefficient  technique 
(Mitalas  and  Stephenson,  1967).  The  fourth  subroutine  is  used  repeatedly 
by  the  first  three.  The  procedure  for  inputting  data  into  the  program 
is  given  in  Appendix  F  and  a  complete  program  listing  is  found  in  Appendix 

This  program  allows  the  user  to  choose  the  time  increment  (At)  at 
which  he  desires  to  run  the  program.  The  only  restriction  is  that  the 
selected  increment  must  divide  evenly  into  24.  The  smallest  time  incre- 
ment, and  the  most  accurate,  would  be  a  one  hour  interval.  The  accuracy 
of  any  interval  greater  than  four  is  suspect. 

The  program  user  also  has  the  option  of  selecting  the  hour  of  the 
day  at  which  to  initiate  load  calculations.  This  feature  is  useful  when 
using  time  increments  greater  than  one  hour.  The  optimal  starting  time 
can  be  chosen  to  produce  cooling  loads  which  will  mostly  closely  approxi- 
mate those  obtained  by  hourly  calculations. 

The  main  program,  a  flow  chart  of  which  is  shown  in  Figure  4,  serves 
to  input  the  majority  of  the  data  and  also  calls  the  three  option-specific 
subroutines.  An  additional  function  of  the  main  program  is  to  calculate 
the  enthalphy  (H)  of  the  outside  air  using  a  mathematical  model  of  the 
psychometric  chart  (Brooker,  1967).  The  formulation  is  divided  into  two 
temperature  regions.  For  temperatures  greater  than  32°F,  H  is  calculated 
as  described  in  equation  9.  And  for  temperatures  less  than  or  equal  to  32°F, 
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Figure  4.  Main  Program 
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Figure  4.      (continued) 
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H  is  calculated  by  equation  10, 

H  -  0.2405  (T-459.69)  +  X  {0.448  T  -  0.1783  T  +  864.7168}  (9) 

H  -  0.2405  (T-459.69)  +  X  {0.448  T  -  0.01377T  +  862.3629}  (10) 

where  T  is  equal   to  the  dry  bulb  temperature  in  degrees  Rankine.     The 
quantity  x  is  determined  by  equation  11, 

X=0-6219       14.696-Pv  (ID 

where  P     is  the  actual   vapor  pressure  of  the  air  at  a  given  temperature 
and  relative  humidity.     Then,  P     is  calculated  by  equation  12. 

P     =  P     •   RH  (12) 

v        s 

In  this  equation  RH  is  the  relative  humidity  in  decimal  units  of  the  out- 
side air  and  P  is  the  saturation  vapor  pressure  of  air  for  a  given  tempera- 
ture (T).  P  may  be  calculated  by  equation  13  or  14,  as  appropriate. 

Ps  =  exp  (54.6329  -  1230j-688  -  5.16923  In  T)   T>32°F    (13) 

Ps  =  exp  (23.3924  -  1128^6489  -  0.46057  In  T)    T<32°F    (14) 

The  main  program  also  prints  a  large  portion  of  the  inputted  climatic 
data  so  that  its  accuracy  can  be  confirmed. 

A  flow  chart  describing  the  subroutine  CLDOPT  is  shown  in  Figure  5.  After 
reading  the  necessary  data,  CLDOPT  repeatedly  calls  TRANEX  until  the  exter- 
ior transmission  load  has  been  calculated  for  each  outside  wall.  CLDOPT  then 
performs  the  same  calculations  for  the  roofs  after  reading  the  appropriate 
new  data.  The  remaining  calculations  are  done  in  a  \zery   straightforward 
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Figure  5.     (continued) 
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manner.     The  units  in  all   intermediate  calculations  are  BTU/day,  choosing 
the  twenty-first  day  of  each  month  as  representative.     The  final  results 
are  given  in  terms  of  Btu  per  pound  of  finished  product. 

Flow  charts  of  the  subroutines  for  hot-processing  options  I  and  II, 
OPT!  and  0PT2  respectively,  have  not  been  made  because  of  their  inherent 
similarity  to  CLDOPT.     Disregarding  the  obvious  changes  in  plant  configura- 
tion,  the  only  difference  between  CLDOPT  and  the  other  option  subroutines 
is  the  inclusion  of  the  shroud  in  the  cold  boning  calculations. 

The  subroutine  TRANEX  is  by  far  the  most  complex  and  its  flow  chart 
is  shown  in  Figure  6.     As  formulated,  the  transfer  function  coefficient 
technique  calculates  the  external   transmission,  q   (t.),   as  function  of  q 

C  I  C 

at  times  prior  to  t.  which  are  initially  unknown.  The  method  used  was  to 
set  all  qg(t)  values  equal  to  zero  at  the  beginning  of  the  calculation. 
An  iterative  procedure  was  employed  which  caused  the  effect  of  this  assump- 
tion to  become  negligible  as  the  calculation  proceeded.  As  programmed, 
TRANEX  calculates  an  exterior  transmission  load  value  for  each  time  incre- 
ment. These  values  are  then  used  to  calculate  a  daily  total.  The  sub- 
routine iterates  on  the  daily  total  values  until  two  adjacent  iterations 
agree  to  within  some  specified  convergence  criterion  which  in  this  case 
was  one  percent.  An  upper  limit  of  20  was  used  on  the  number  of  iterations. 
If  this  limit  is  reached,  an  appropriate  message  is  printed  and  execution 
halted. 
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PLANT  CONFIGURATIONS  AND  SYSTEM  VARIABLES 

The  flow  sequences  of  the  two  hot-processing  options  considered  in 
this  study  are  shown  in  Figures  7  and  8.     These  figures  also  show  the 
associated  mass,  specific  heat  and  entering  and  exiting  temperatures  of  the 
product  contained  in  each  of  the  refrigerated  coolers.     Figure  9  displays 
this  same  information  for  the  conventional   cold-processing  facility  considered 
in  this  analysis. 

The  cooler  dimensions  and  capacities  considered  in  this  study  for 
hot  processing  I,  hot  processing  II  and  cold  processing  are  shown  in 
Figures  10,  11  and  12,  respectively. 

A  listing  of  all  data  requirements  and  the  corresponding  values 
selected  in  this  study  is  given  in  Table  4. 
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RESULTS  AND  DISCUSSION 


A  transient  heat  analysis  was  performed  to  determine  the  rate  of 
product  cooling  under  hot-processing  options   I  and  II  and  cold  processing. 
This  analysis  showed  that  104  hours  was  required  for  the  internal  product 
temperature  to  reach  28.3°F  with  hot-processing  option   I.     It  required  96 
hours  for  the  internal  product  temperature  to  reach  28.3°F  under  hot- 
processing  option  II  and  72  hours  with  conventional   cold-processing 
techniques. 

The  results  from  the  transient  heat  analysis  were  used  in  the  pre- 
viously described  computer  program  to  compute  refrigeration  energy  loads. 
The  computer  program  calculated  the  average  energy  required  for  product 
and  space  refrigeration  for  hot-processing  options  I  and  II  and  cold 
processing. 

A  monthly  comparison  of  the  refrigeration  energy  expended  for  hot- 
processing  options   I  and  II  and  cold  processing  is  shown  in  Figure  13.     A 
comparison  shows  that  in  December,   the  coldest  month,   133.5  Btu  were 
required  to  process  one  pound  (3.1  x  105  Joules  per  kg)  of  beef  by  the 
conventional   cold-processing  technique.     In  the  same  month  93.2  But/lb 
(2.2  x  105  J/kg)  were  required  using  hot-processing  option  I  and  76.9 
Btu/lb  (1.8  x  10s  J/kg)  using  hot-processing  option   II.     This  constitutes 
a  30  percent  savings  for  hot  processing  I,  and  a  42  percent  savings  for 
hot  processing  II  when  compared  to  cold  processing. 

July  was  the  warmest  month  and  results  obtained  for  this  period 
show  that  the  cold-processing  energy  requirements  were  153.8  Btu/lb 
(3.6  x  105  J/kg);  the  hot-processing  I  requirements  were  103.8  Btu/lb 
(2.4  x  105  J/kg);  and  the  hot-processing  II  requirements  were  89.6  Btu/lb 
(2.1  x  10s  J/kg).     The  associated  percentage  energy  savings  in  July  for 
hot  processing  I  was  33  percent,  and  42  percent  for  hot  processing  II. 
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Figure  13.  Refrigeration  Energy  Utilized  on  the 
Twentv-First  Dav  of  Each  Month. 
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Figure  14  shows  the  proportional   amount  of  each  of  the  individual 
system  loads  to  the  total  load  for  each  processing  scheme.     A  comparison 
on  this  basis  emphasizes  the  major  areas  of  energy  savings  between  the  hot- 
processing  options  and  cold  processing.     A  major  contributor  to  the 
hot-processing  savings  was  the  absence  of  the  shroud  load,  inherent  only 
in  cold  processing.     The  elimination  of  the  shroud  alone  generated  30 
percent  of  the  energy  savings  in  hot  processing  I  and  24  percent  of 
the  savings  of  hot-processing  option   II.     A  second  principal   source 
of  energy  reduction  was  the  reduced  product  loads  of  the  hot  options 
caused  by  not  cooling  portions  of  the  carcasses  that  will  be  rendered. 
The  reduced  product  load  constituted  29  percent  of  the  energy  reduc- 
tion between  hot  processing  I  and  cold  processing  and  32  percent  of 
the  reduction  between  hot  processing  II  and  cold  processing.     External 
transmission  was  another  significant  area  of  energy  savings.     It  was 
responsible  for  21  percent  of  the  energy  savings  associated  with  hot 
processing  I  and  9  percent  of  the  hot-processing  II  reduction.     The 
reduced  electrical   and  lighting  loads  of  hot  processing  II  were  signi- 
ficant sources  of  energy  savings  for  this  option.     These  savings  are 
the  result  of  the  greatly  reduced  cooler  volume  associated  with  hot 
processing  II. 

Several   sample  calculations  were  performed  to  determine  the 
effect  of  changing  the  time  interval  between  calculations  from  one 
to  four  hours.     The  effect  of  starting  time,  when  using  four  hour 
intervals,  was  also  investigated.     The  change  in  the  external   trans- 
mission load,  using  any  of  these  procedures,  was  less  than  2  percent, 
thus,  not  practically  significant. 
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SUMMARY  AND  CONCLUSIONS 

A  transient  heat  analysis  indicated  that  104  hours  was  required 
for  the  internal   product  temperature  to  reach  28.3°F  (-2.1°C)  when  the 
carcass  was  conditioned  for  8  hours  at  60°F  (15.6°C)  and  subsequently  hot 
processed  into  subprimals,  vacuum  packaged  and  boxed.     The  analysis  also 
showed  that  96  hours  was  required  for  the  internal  product  temperature 
to  reach  28.3°F  (-2.1 °C)  when  the  carcass  was  fabricated  into  subprimals 
and  vacuum  packaged  and  boxed  immediately  following  slaughter.     Conven- 
tionally processed  beef  required  72  hours  to  reach  an  internal   tempera- 
ture of  28.3°F  (-2.1°C). 

A  computer  program  was  developed  to  calculate  the  average  refriger- 
ation energy  requirements  of  both  hot-processing  schemes  and  cold  process- 
ing.     A  listing  is  shown  in  Appendix  G. 

An  energy  comparison  among  the  two  hot-processing  option  and  cold 
processing  indicated  that  a  significant  energy  savings  may  be  achieved 
by  utilizing  hot  processing  in  the  beef  industry.     This  requires  adjust- 
ing existing  plants'    refrigeration  and  holding  facilities  or  constructing 
new  facilities  to  accomodate  hot  processing. 

Beef  carcasses  conditioned  at  60°F  (15.6°C)   for  eight  hours,  and 
immediately  fabricated  into  subprimals,   vacuum  packaged  and  boxed  required 
approximately  50  Btu  per  lb  of  finished  product  (1.2  x  10s  Joules  per  kg) 
less  energy  than  carcasses  that  were  conventional   processed  into  subprimals 
after  chilling  for  72  hours.     This  is  a  32  percent  reduction  in  the  total 
refrigeration  energy  expended  during  product  fabrication. 

Carcasses  which  were  fabricated  into  subprimals,  and  vacuum  packaged 
and  boxed  immediately  following  slaughter  required  approximately  64  Btu 
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per  lb  of  finished  product  (1.5  x  10     Joules  per  kg)  less  energy  than 
cold-processing  carcasses.     This  is  a  42  percent  reduction  in  the  refrig- 
eration energy  required  under  the  cold-processing  technique. 

These  energy  comparisons  were  based  upon  model  facilities  that 
slaughter  sixty  head  of  cattle  an  hour  for  eight  hours  a  day.     However, 
conclusions  found  for  this  size  hot-processing  facility  should  not  be 
unreservedly  ascribed  to  other  sized  plants. 
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DIRECTIONS  FOR  FUTURE  RESEARCH 

Hot  processing  is  an  effort  to  minimize  the  product  load  component 
of  the  total   refrigeration  energy  requirement.     Further  research  could 
extend  load  reduction  and  optimization  to  other  components  of  the  refrig- 
eration load.     A  prime  consideration  for  future  efforts  would  be  to 
determine  the  minimum  amount  of  cooler  space  required  to  efficiency 
process  a  beef. 

The  energy  savings  due  to  hot  processing  is  only  one  of  the  poten- 
tial economic  advantages.     An  overall   economic  analysis  to  include  labor, 
capital   investment,  etc.   is  required  to  evaluate  the  full   potential  of 
hot  processing.     A  cost-benefit  analysis  could  weigh  the  economic  bene- 
fits of  hot  processing  against  any  detriments.     These  analyses  could 
further  aid  the  industry  in  making  the  decision  to  implement  hot  process- 
ing. 
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APPENDIX  A 
TRANSIENT  HEAT  ANALYSIS  OF  HOT-PROCESSING  OPTION  I 


Problem:  To  evaluate  the  internal  carcass  temperature  and  the  average 

carcass  temperature  after  8  hours  in  the  60°F  conditioning  room. 

Assumptions:  Consider  the  beef  round  to  be  represented  by  the  shape  of 
a  rectangle.  It's  dimensions  are  20"  x  14.5"  x  8.5" 

Solution:  Perform  a  transient  heat  analysis  on  the  beef  round  using 
Heisler's  (1947)  equation  15  and  chart  (Figure  15). 


Total 
Round 


T-To 


T  -TQ 


14.5" 
slab 


T  =  final  internal  temperature,  F 
T  =  uniform  initial  temperature,  F 
Too  ■  surrounding  fluid  temperature,  F 


a8 
To  use  Heisler's  chart  must  know  t-j-  and  K/hL  where 

L  ■  1/2  slab  thickness,  ft 

9  ■  time  elasped,  hours 

a  =  thermal  diffusivity,  sq  ft  per  hour 

K  =  thermal  conductivity,  Btu  per  (hr)  (ft)  (deg  F) 

h  =  surface  heat  transfer  coefficient,  Btu  per  (hr)  (sq  ft)  (deg  F) 

a  =  0.0048  sq  ft  per  hour  (ASHRAE  Fundamentals,  1977d) 

h  =  2.05  Btu  per  (hr)  (sq  ft)  (deg  F)  (ASHRAE  Fundamentals,  1977d) 


*t* 


eat 


0.28  Btu  per  (hr)   (ft)   (deg  F)   (ASHRAE  Fundamentals,  1977d, 

Lentz,  1961) 
From  Fig.    15,  for  0=8  hrs. 
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Figure  15.     Miripiane  temperature  for  an  infinite  plate 
of  thickness  2L,  from  Heisler  (1947). 
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After  8  hours  at  60°F,  with  an  initial  temperature  of  105°F  (ASHRAE  Appli- 


cations, 1974a). 


T-60 


105-60 


=  (1)  (0.99)  (0.80) 


T  =  95.6UF 


To  determine  the  average  temperature  of  the  beef  round  after  8  hours  in 
60°F  room  use  Heisler's  chart,  Figure  16  which  relates  temperature 
at  any  position  as  a  function  of  center  temperature. 


slab 


K  /hL 


X/L 


T(x,0rT" 
T(o,9)-T» 


20" 

0.164 

0.5 

0.78 

14.5" 

0.226 

0.5 

0.80 

8.5" 

0.386 

0.5 

0.85 

Now  using  equation 

16 

(Heisler,  1947) 

(Tx,0-O 
lTo..-T-J 

Total 

., 

To,9'T=o  [  20" 
L      J  sla 

J 

X 

T0 

,0 

,8 

■T 

00 

-T-J 

14.5" 
slab 

[Tx,o"T. 
iT0,9-T- 

T    =  average  temperature  of  slab  for  X/L  =  0.5,  T 

A  jU 

TM   =  surrounding  air  temperature,  F 

T  Q  ■  center  temperature,  F  (previously  determined) 

O,o 


(16) 
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After  8  hours  in  50  F  room 
■60 
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95.6-60 
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T. 


x,o 


78.3°F 
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Figure  16.  Temperature  as  a  Function  of  Center  Temperature  in  an 
Infinite  Plate  of  Thickness  2L,  from  Heisler  (1947). 
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Problem: 


Solution: 


To  evaluate  the  time  required  for  an  essentially  solid  box 
(9.5"  x  22"  x  17")  of  meat  product  to  reach  an  internal 
temperature  of  28.84 of  in  the  28°F  shipping  cooler. 

Perform  a  transient  heat  analysis  using  Heisler's  (1947) 
equation  15  and  chart  (Figure  15).   Must  reiterate,  varying 
time  (e)  until  desired  temperature  is  reached. 


T-To 


Total 
box 


o  °° 


22" 
slab 


T-To 

IV1 


17" 
slab 


[; 


T-To, 


(15) 


9.5" 
slab 


o 

Too 


final  internal  temperature,  deg  F 
initial  uniform  temperature,  deg  F 

surrounding  fluid  temperature,  deg  F 


To  use  Heisler's  chart  must  know  a6/L2  and  K/hL  where 

L  =  1/2  slab  thickness,  ft 

8  =  time  elasped  in  cooler,  hours 

a  =  thermal  diffusivity,  sq  ft  per  hour 

K  =  thermal  conductivity,  Btu  per  (hr)  (ft)  (deg  F) 

h  =  surface  heat  transfer  coefficient,  Btu  per  (hr)  (sq  ft)  (deg  F) 

a  =  0.0048  sq  ft  per  hour(ASHRAE  Fundamentals,  1977d) 

h  =  2.05  Btu  per  (hr)  (sq  ft)  (deg  F)  (ASHRAE  Fundamentals,  1977d) 

K    =  0.28  Btu  Der  (hr)  (ft)  (deg  F)  (ASHRAE  Fundamentals,  1977d  and 
meat 


:< 


box 


Lentz,  1961) 

1.74  Btu  per  (hr)  (ft)  (deg  F) 


weight  box  =  3.5  pounds,  weight  meat  =  90  pounds 


K 


90  lb 


Btu 


composite 


Iotf      °-28  wmm 


From  Figure  15;  for  e  =  96  hours 
slab  ae/L2 


K/hL 


L.74  Btu       I 
(hr)(ft)(Fj  f 


T-Tcq 


0.33  Btu 
(hr)(ft)(F) 


22" 
17" 
9.5" 


0.548 
0.920 
2.940 


0.176 
0.227 
0.407 


0.55 
0.30 
0.027 
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After  96  hours  in  a  28°F  room,  with  an  initial  temperature  of  95.6  F 


1^28 
95 


§328  ■  (°-55)  (°-3°)  (°-027) 


T  =  28.30  *  28.84UF 

To  determine  the  average  temperature  of  the  box  after  96  hours  in  the 
28°F  room  use  Heisler's  (1947)  chart,  Figure  16,  which  relates  the 
temperature  at  any  position  as  a  function  of  center  temperature. 


slab 


22" 

17" 
9.5" 


K  /hi 


X/L 


0.176 

0.5 

0.227 

0.5 

0.407 

0.5 

T(x,orT°° 
T(0,9)"T- 


0.75 
0.83 
0.85 


Now  using  equation  16  (Heisler,  1947) 


x,o"  °° 
o,e~  oo 


=  (Ixifii-i     (!** 


-T 


total 


To,e'Tcof  22"   I To,e~Too  j  17" 


slab 


slab 


Tx,o"T 

T0,9-T»J 

9.5" 
slab 

(16) 


where: 
x,o 


0,9 


average  temperature  of  slab  for  X/L  =0.5,  deg  F 

surrounding  air  temperature,  deg  F 

center  temperature,  deg  F  (previously  determined) 


Tx,o"28   =  (0.75)  (0.83)  (0.85) 
28.30-28 


TXj0=28.16"F 
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Problem: 


Solution: 


APPENDIX  B 
TRANSIENT  HEAT  ANALYSIS  OF  HOT-PROCESSING  OPTION  II 


To  evaluate  the  time  required  for  an  essentially  solid  box 
(9.5"  x  22"  x  17")  of  packaged  meat  product  to  reach  an 
internal  temperature  of  28.84 °F  in  the  28°F  shipping  cooler. 

Perform  a  transient  heat  analysis  using  Heisler's  (1947) 
equation  15  and  chart  (Figure  15). 


Total 
box 


22" 
slab 


J  T-Too  1        I  T-Tco  I 

\ryrr  „„      \r^r[ 


17" 
slab 


9.5" 
slab 


T   =  final  internal  temperature,  deg  F 
T   =  initial  uniform  temperature,  deg  F 
Ta>  =  surrounding  fluid  temperature,  deg  F 


To  use  Heisler's  chart  must  know  rf-  and  j-p  where 


L  =  1/2  slab  thickness,  ft 

9  =  time  elasped  in  cooler,  hours 

a  =  thermal  diffusivity,  sq  ft  per  hour 

K  =  thermal  conductivity,  Btu  per  (hr)  (ft)  (deg  F) 

h  =  surface  heat  transfer  coefficient,  Btu  per  (hr)  (sq  ft)  (deg  F) 


(15) 


a   =  0.0048  sq  ft  per  hour  (ASHRAE  Fundamentals,  1977d) 
h   =  2.05  Btu  per  (hr)  (sq  ft)  (deg  F)  (ASHRAE  Fundamentals,  1977d) 
Kmeat  =  0.28  Btu  per  (hr)  (ft)  (deg  F)  (ASHRAE  Fundamentals,  1977d  and 
Lentz,  1961) 

Kbox  =  L74  Btu  per  (hr)  (ft)  (deg  F) 


Weight  box  =  3.5  pounds,  weight  meat  =  90  pounds 


composite 


90  lb 
93.5  lb 


1 


0.28 


Btu 


IhrUftKF) 


3.5  lb 
93.5  lb 


1.74  Btu 
(hr)(ft)(Fj 
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0.33  Btu 
THF)(ft)(F) 


From  Figure  15;  for  9  =  96  hours 
slab  a6/L 


T-T« 


K/hL 


VT~ 


0.55 
0.30 
0.027 

After  96  hours  in  a  28°F  room,  with  an  initial   temperature  of  105  F 


22" 

0.548 

0.176 

17" 

0.920 

0.227 

9.5" 

2.940 

0.407 

T-28 
105-28 


(0.55)    (0.30)    (0.027) 


T  =  28.343  ■  28.84UF 


To  determine  the  average  temperature  of  the  box  after  96  hours  in  the 
28°F  room  use  Heisler's  (1947)  chart,  Figure  16,  which  relates  the 
temperature  at  any  position  as  a  function  of  center  temperature. 


slab 

K  /h  L 

X/L 

T(x,orT» 
T(o,e)-T» 

22" 
17" 
9.5" 

0.176 
0.227 
0.407 

0.5 
0.5 
0.5 

0.75 
0.83 
0.85 

Now  using  equation  16.  (Heisler,  1947) 


where: 
x,o 

CO 

To,e  = 


total 


average  temperature  of  slab  for  X/L  =  .5,°F 

surrounding  air  temperature,  °F 

center  temperature,  °F  (previously  cetejyiined) 


x,o~ 
To,e"T« 


(16) 


9.5" 
slab 


After  96  hours  in  28  F  room 


T   -28 
x,o 

28.34-28 


=  (0.75)  (0.83)  (0.85) 


Tx   =  28.18UF 
x,o 
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APPENDIX  C 
TRANSIENT  HEAT  ANALYSIS  OF  COLD  PROCESSING 


Problem:  To  evaluate  the  internal  carcass  temperature  and  the  average 
carcass  temperature  after  24  hours  in  the  30  F  chill  cooler. 

Assumptions:  Consider  the  beef  round  to  be  represented  by  the  shape 
of  a  rectangle.  Its  dimensions  are  20"  x  14.5"  x  8.5" 


Solution:  Perform  a  transient  heat  analysis  on  the  beef  round  using 
Heisler's  (1947)  equation  (15)  and  chart  (Figure  15  ). 

T-Too         _  I  T-Too  [       I  T-T< 

20" 
slab 


VT«   Total      VT~  f  20"    1  VTa 


Round 


T  =  final  internal  temperature,  F 
T  =  uniform  initial  temperature,  F 
^  ■  surrounding  fluid  temperature,  F 


To  use  Heisler's  chart  must  know  ff-  and  K/hL  where 

L  =  1/2  slab  thickness,  ft 

6  =  time  elasped,  hours 

a  =  thermal  diffusivity,  sq  ft  per  hour 

K  =  thermal  conductivity,  Btu  per  (hr)  (ft)  (deg  F) 

h  =  surface  heat  transfer  coefficient,  Btu  per  (hr)  (sq  ft)  (deg  F) 

a  =  0.0048  sq  ft  per  hour  (ASHRAE  Fundamentals,  1977d) 
h  =  2.05  Btu  per  (hr)  (sq  ft)  (deg  F)  (ASHRAE  Fundamentals,  1977d) 
Kmeat  ■  0.28  Btu  per  (hr)  (ft)  (deg  F)  (ASHRAE  Fundamentals,  1977d  and 
Lentz,  1961) 


From  Fig.  15,  for  9  = 


24  hours 


59 


T-T 


b  IdD 

axil  U 

iv  ni_ 

1>T 

O       00 

20" 

0.166 

0.164 

0.97 

14.5" 

0.316 

0.226 

0.75 

8.5" 

0.918 

0.386 

0.38 

After  24  hours  in  30°F  room,  with  an  initial  temperature  of  105°F 

]^^    =  (0.97)  (0.75)  (0.38) 

T  =  50.73°F 

To  determine  the  average  temperature  of  the  round  after  24  hours  in  the 
30°F  chill  cooler  use  Heisler's  (1947)  Chart,  Figure  16,  which  relates 
temperature  at  any  position  as  function  of  center  temperature. 


slab 


K/hL 


X/L 


T(x,orTc 

T(o,9)-Tc 


20" 

14.5" 
8.5" 


0.164 
0.226 
0.386 


0.5 
0.5 
0.5 


0.76 
0.80 

0.85 


Now  using  equation  16  (Heisler,  1947) 


°'9  °°  i  Total 


T   -T 

X,0   oo 

To,e"To=|  20" 
slab 


T   -T 
x,o  °° 

To,9"T=o  I  14.5" 
slab 

_o 


T    =  average  temperature  of  slab  for  X/L  =  0.5,  F 

A  )U 

T^   =  surrounding  air  temperature, °F 

T  Q  =  center  tempera ture,°F  (previously  determined) 
o  ,y 

After  24  hours  in  30°F  room 

Tx  o"30 

5-g;73_30  »  (0.76)  (0.80)  (0.85) 


(w 


T0,8-To 


(16) 


8.5" 
slab 


Tx,o  =  40-71°F 
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Problem:  To  evaluate  the  internal  carcass  temperature  and  average  carcass 
temperature  after  24  and  48  hours  in  the  32°F  holding  cooler. 

Assumptions:  Consider  the  beef  round  to  be  represented  by  the  shape  of 
a  rectangle.  Its  dimensions  are  20"  x  14.5"  x  8.5". 

Solution:  Perform  a  transient  heat  analysis  on  the  internal  round  using 
Heisler's  (1947)  equation  15  and  chart  (Figure  15  ). 


fet} 


T-T0 


T-Ta 


Total 
Round 


T  -T 

O   =0 


20" 
slab 


T  =  final  internal  temperature,  °F 
TQ  =  uniform  initial  temperature,  °F 
^  =  surrounding  fluid  temperature, °F 


14.5" 
slab 


k 


vt» 


(15) 


8.5" 
slab 


To  use  Heisler's  chart  must  know  £§■  and  K/hL  where 

L  =  1/2  slab  thickness,  ft 

6  =  time  elasped,  hours 

a  =  thermal  diffusivity,  sq  ft  per  hour 

K  =  thermal  conductivity,  Btu  per  (hr)  (ft)  (deg  F) 

h  =  surface  heat  transfer  coefficient,  Btu  per  (hr)  (sq  ft)  (deg  F) 

a  =  0.0048  sq  ft  per  hour  (ASHRAE  Fundamentals,  1977d) 

h  =  2.05  Btu  per  (hr)  (sq  ft)  (deg  F)  (ASHRAE  Fundamentals,  1977d) 

Vat  =  °'28  Btu  per  (hr)  (ft)  (de9  F)  (ASHRAE  Fundamentals,  1977d  and 
Lentz,  1961) 

From  Fig.  15,  for  6  ■  24  hours. 

slab        ae/L2        K/hL 


r-T 


I  -  loo 

VTco 
0 


20" 
14.5" 
3.5" 


0.166 
0.316 
0.918 


0.164 
0.226 
0.386 


0.97 
0.75 
0.38 


After  24  hours  in  32°F  holding  cooler,  with  an  initial  temperature  of  50.73°F. 


T-32 


50.;3.32  =  (0-97)  (0'75)  (0'38) 


37.17°F 
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After  48  hours  in  32°F  holding  cooler,  with  an  initial  temperature  of  50.7°F 

l~732l7_Z2    =  (°-97)  (°-75)  (°'38) 

T  =  33.43°F 

To  determine  the  average  temperature  of  the  round  after  24  and  48  hours 
in  the  32°F  holding  cooler  use  Heisler's  (1947)  chart,  Figure  ^,   which 
relates  temperature  at  any  position  as  a  function  of  center  temperature. 


slab 

K/hL 

X/L 

T(x,orT~ 
T(o,9)-T~ 

20" 
14.5" 
8.5" 

0.164 
0.226 
0.386 

0.5 
0.5 
0.5 

0.76 
0.80 
0.85 

Now  using  equation  16  (Heisler,  1947) 

Tx,o"T°°  I  _  J  Tx,o"T°°  I  J  x,o"  ° 

T  ft-T  f  ,  t  ,      To.9"T°°l  20"  lTo,6"T«>   14.5" 

0.8  -J  Total     L      J  slab  L      J  slab 


Or 


T    =  average  temperature  of  slab  for  X/L  =  0.5,UF 

T    =  surrounding  air  temperature , T 

T  fl  ■  center  temperature ,°F  (previously  determined) 

0,9 

After  24  hours  in  32°F  holding  cooler 

T   -32 

3^;°7_32  =  (0.76)  (0.80)  (0.85) 

T    =  34.67°F 
x,o 

After  48  hours  in  32°F  holding  cooler 


T   -32 
x,o 

33.43-32 

=  (0.76)  (0.30)  (0.85) 

T    =  32.74°F 
x, 
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Problem:  To  evaluate  the  internal  and  average  temperature  of  an  essen- 
tially solid  box  (9.5"  x  22"  x  17")  of  packaged  meat  product 
after  24  and  48  hours  in  the  28°F  shipping  cooler. 

Solution:  Perform  a  transient  heat  analysis  using  Heisler's  (1947) 
equation  15  and  chart  (Figure  15  ). 


total 
box 


T-Too 


T-T0 


22" 

slab 


o  °° 


17" 
slab 


T-Ta 


VT° 


9.5" 
slab 


(15) 


T  =  final  internal  temperature,  deg  F 
T  =  initial  uniform  temperature,  deg  F 

Too  =  surrounding  fluid  temperature,  deg  F 


To  use  Heisler's  chart  must  know  a9/L2  and  K/hL  where 

L  =  1/2  slab  thickness,  ft 

9  »  time  elasped  in  cooler,  hours 

a  =  thermal  diffusivity,  sq  ft  per  hour 

K  =  thermal  conductivity,  Btu  per  (hr)  (ft)  (deg  F) 

h  =  surface  heat  transfer  coefficient,  Btu  per  (hr)  (sq  ft)  (deg  F) 

a  =  0.0048  sq  ft  perhour  (ASHRAE  Fundamentals,  1977d) 

h  =  2.05  Btu  per  (hr)  (sq  ft)  (deg  F)  (ASHRAE  Fundamentals,  1977d) 

Kmeat  =  0-28  Btu  per  ^hr^  ^ft^  ^deg  F^  (ASHRAE  Fundamentals,  1977d  and 
Lentz,  1961) 


"box 


=  1.74  Btu  per  (hr)  (ft)  (deg  F) 


weight  box  =  3. 

5  pound 

s,  weight  meat  =  90  pounds 

K                =  1  90 
composite      ]    93 

.5  lb  f 

f°-2S  AtFI  j  *  j 

93.5 

lb 
lb 

jl.74  Btu 

f      i(iir)(ft)(F) 

From  Figure  15 
slab 

a  /L                      K/hL 

T-Too 

O          00 

22" 

0.137                    0.176 

0.99 

17" 

0.230                    0.227 

0.88 

9.5" 

0.735                    0.407 

0.45 

0.33  Btu 
(hr)(ft)(F) 
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After  24  hours  in  a  28°F  room,  with  an  initial  temperature  of  33.43  F 


T-28 


33.43-28 


=  (0.99)  (0.88)  (0.45) 


T  =  30.13UF 


After  48  hours  in  28  F  room 


T-28 

30. 


^g  =  (0.99)  (0.88)  (0.45) 


T  =  28.84UF 

To  determine  the  average  temperature  of  the  box  after  24  and  48  hours  in 
the  28°F  room  use  Heisler's  (1947)  chart,  Figure  16,  which  relates  the 
temperature  at  any  position  as  a  function  of  center  temperature. 


slab 

K/hL 

X/L 

• 

T(x,o)_T- 
T(o,erT~ 

22" 
17" 
9.5" 

0.176 
0.227 
0.407 

0.5 

0.5" 

0.5 

0.75 
0.83 
0.85 

Now  using  equation  16  (Heisler,  1947) 


Tx,o-T° 
lT0,9"T°°- 


Total 


x,o"  OT 
To,9-To 


20" 
slab 


T    =  average  temperature  of  slab  for  X/L  =  0.5,  F 
T    =  surrounding  air  temperature,  F 

00 

T    =  center  temperature,  F  (previously  determined) 

0,6 


After  24  hours  in  28  F  room 

(0.75)  (0.83)  (0.85) 


Tx,o-28 
30.13-28 


x,o 


29.18°F 


After  48  hours  in  28  F  room 


T   -28 

2|^g  ■  (0.75)  (0.83)  (0.85) 

T    =  28.46°F 
x,o 
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APPENDIX  D 
SAMPLE  CALCULATION 

External  Transmission   East  wall  on  July  21,  1976 


<e.r  =  A 


Z  b-  (ta  m   nA)  "  s  dn  VMe,r-nA'  "  t 


jn=o 


(qe,r-nA)  -  t   E  c  I      (4) 


n  '"e.r-nA'   n=1  n  jj re  n=Q 


where:  qQ  _  ■  heat  gained  by  the  space  in  Btu  per  hour 
e  ,r 

A  =  exposed  area,  sq  ft 

A  ■  time  interval  in  hours 

t  =  sol -air  temperature,  deg  F 

t  ■  constant  indoor  air  temperature,  deg  F 

b  ,c  ,d  =  transfer  function  coefficient  characteristic  of  wall 
n  n  n 

construction 

For  this  example: 

A   =  92  ft  x  28.5  ft  =  2,622  ft 

trc  =30°F 

b  ,c  ,d  =  transfer  functions  for  2  in.  insulation  on  8  in. 

l.w.  concrete  block  from  ASHRAE  Fundamentals  (1977) 
n=o     n=l      n=2       n=3      n=4      n=5 
b   0.0001   0.00081   0.00339    0.00196   0.00018 
d   1.0000  -1.52490   0.69509   -0.11032   0.00399   -0.00001 

Z  cn   =  0.00635 

<*I*         CAD 

0        0 

where:  t   =  outdoor  air  temperature,  deg  F 
a   =  absorbance  of  surface 
I.  =  total  solar  radiation  incident  upon  surface,  Btu  per  (hour) 

(sq  ft) 
h   =  coefficient  of  heat  transfer,  Btu  per  (hour)  (sq  ft)  (deg  F) 
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E   =  hemispherical  emittance  of  surface 

AR  =  difference  between  longwave  radiation  incident  upon  the 
surface  from  the  sky  and  surroundings  and  the  radiation 
emitted  by  blackbody  at  outdoor  air  temperature,  Btu  per 
(hour)  (sq  ft) 


For  this  example: 


a/h  =0.30  for  dark-colored  surface  from  ASHRAE  Fundamentals  (1977) 

o 
AR   =0  for  vertical  surface  from  ASHRAE  Fundamentals  (1977) 

E   =  0.82  for  masonry  materials  from  ASHRAE  Fundamentals  (1977a) 
h   =  4.0  for  7.5  mph  wind  from  ASHRAE  Fundamentals  (1977a) 
It   =  1.15  (SHGF) 
SH6F  is  a  function  of  latitude,  date,  and  solar  time 


Fc 

r  July  21  at 
Solar  Time 

40  deq,  N 

SHGF- 
(Btuh  pe 

orth 

Latitude,  East  ! 

Exposure 

CST-DST 

/ 
r  ft: 

,,     !t      «*«•" 

W°F> 

8:00  A.M. 

6:30  A.M. 

170.5 

196.06     76 

134.8 

12:00  A.M. 

10:30  A.M. 

113.5 

130.53     85 

124.2 

4:00  P.M. 

2:30  P.M. 

169.5 

194.93     87 

145.5 

8:00  P.M. 

6:30  P.M. 

69.5 

79.93     75 

99.0 

12:00  P.M. 

10:30  P.M. 

0 

0     72 

72 

4:00  A.M. 

2:30  A.M. 

0 

0     78 

78 

Therefore: 

—                 — 

<L  „  -  2622  < 
e,r 

f 

0.0001(134.8) 
+0.00081(78) 
+0.00339(72) 
+0.00196(99) 
+0.00018(145.5) 

- 

1.0000(0) 
-1.52490(0) 
+0.69509(0) 
-0.11032(0) 
+0.0039  (0) 
-0.00001(0) 
- 

- 

30(0.00635) 
J 

qa   =  918.9  Btu  per  hr 
e  ,r 


Repeat  this  calculation  in  successive  cycles  to  nullify  effect  of  initially 


setting  q    ,  equal  to  zero 
e,r-nA 


1/  From  ASHRAE  Fundamentals  (1977b) 
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Interior  Transmission 

q  =  U  A  (tb  -  trc)  (6) 

where:  q  =  heat  gained  by  the  space  in  Btu  per  hour 

U  =  overall  coefficient  of  heat  transfer,  Btu  per  hr  (sq  ft)  (deg  F) 

A  ■  exposed  area,  sq  ft 

tb  =  temperature  of  the  adjacent  space,  deg  F 

t   =  temperature  of  the  cooler  space 

For  this  example: 

A  =  48  ft  x  28.5  ft  =  1,368  ft2 

V  70°F 

trc  -  30°F 

U  =  0.0919  Btu  per  hr  (sq  ft)  (deg  F)  for  wall  constructed  with 
4  inl.w.  concrete  block  with  0.75  in  plaster 

Therefore:  q  ■  0.0919  hr^  ^  (1,368  ft2)  (70°F  -  30°F) 
q  =  5028.8  jjjjH. 

Infiltration 

q  =  4.5  (cfm)  (Ah)  (7) 

where:  q  =  quantity  of  heat  load,  Btu  per  hour 

cfm  =  volume  of  air  entering  the  space,  cubic  feet  per  minute 
Ah  =  difference  in  enthalphy  between  inside  and  outside  air, 
Btu  per  pound  dry  air 

For  this  example  using  July  mean  temperatures  and  relative  humidity: 

houtside  air  '■   38"°  Midni9ht  CST 

K~si—  a     9-6  For  28°F  at  85%  R.H. 

cooler 

cfm  =  14  cu  ft  per  min  from  ASHRAE  Fundamentals  (1977c) 
Therefore: 

"infiltration4-5'14)  (38-9.6)  •  1.789.2  fSi 
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Personnel 

q  ■  N  (Q) 

where:  q  =  heat  gained  by  the  space  in  3tu  per  hour 

N  =  number  of  persons  entering  refrigerated  space 

Q  =  heat  equivalency  per  person,  Btu  per  hour 

For  this  example: 

K  •  7  persons 

Q  =  950  Btu  per  hr  (in  32°F  cooler) 


Therefore: 


q  =  7  person  (950  ^  ) 
q  "  6,650^- 


Product  Load   480  carcasses 

q  =  mc  (tx  -  t2) 

where:  q  =  quantity  of  heat  to  be  removed,  Btu 

m  =  weight  of  product,  lbs 

c  ■  specific  heat  of  carcass,  Btu  per   (lb)  (deg  F) 

t.  ■  temperature  carcass  as  it  enters  cooler 

t2  =  temperature  carcass  as  it  exits  cooler 

For  this  example: 

m  =  480  carcass  (520  lb  per  carcass)  =  297,500  lb 

c  ■  0.75  Btu  per  (lb)  (deg  F) 

t,  ■  98.5°F 

t2  =  A0.7°F 

Therefore: 

q  =  197,500  lb  (0.75  |^r)  (98.5  -  43.7°F) 

q  =   3.556  x  107  Btu 


7^ 


Shroud  Load 

qtotal 

^sensible   qlat 

^sensible 

=  mn  cw  (W 

qlatent 

hfg  (mw} 
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(2) 

(3) 

where:  qf  .  ,  =  total  heat  gained  by  sensible  and  latent  factions,  Btu 

m  =  total  weight  of  water  loss,  lb 

c  =  specific  heat  of  water,  Btu  per  (lb)  (deg  F) 

t  =  temperature  of  shroud  water,  F 

t  =  temperature  of  chill  cooler,  F 

h-  ■  heat  of  vaporization  of  water,  Btu  per  lb  water 

For  this  example: 

m   =  6  pounds  water  per  carcass  (480  carcasses)  =  2,880  lb 

cw  =  1.0  Btu  per  (lb)  deg  F) 

t   =  125°F 

tr  =  30°F 

t   =  1049.5  Btu  perlb  water 

Therefore: 

Sensible     =     273,600  Btu 

latent    =     3,022,560  Btu 

qtotal  =  3>296>160  Btu 

Electrical  -  Motors 

For  24  horsepower  motor  with  motor  losses  outside  the  refrigerated  space, 
q  =  24  hp  (2,545  jj^)  =  61,080  ^ 

2,545  conversion  given  in  ASHRAE  Fundamentals  (1977c) 
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Electrical  -  Lighting 
q  =  total  wattage  factor  x  use  factor  x  special  allowance  factor  x  3.41   (8) 

For  this  example: 

total  wattage  =  3  watt/ft2  (8,832  ft2)  =  26,496  watt 

use  factor  =  1 

special  allowance  factor  =  1 

q  =  26,496  watt  (1)  (1)  3-w^tu/hr   =  90351.4  Btu/hr 
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APPENDIX  E 
SOLAR  TIME 

Actual  solar  time  may  be  calculated  using  equation  15  from  ASHRAE 
Applications,  1974b. 

AST  =  Local  Standard  Time  +  Equation  of  Time  +  4  x  (number  of    (15) 
minutes  east)  or  -  4  x  (number  of  minutes  west)  of  the 
Local  Standard  Time  Meridian 

The  longitudes  of  the  six  Standard  Time  Meridians  which  affect  the  U.S. 

are:  Eastern  ST,  75  deg;  Central  ST,  90  deg;  Mountain  ST,  105  deg; 

Pacific  ST,  120  deg;  Yukon  ST,  135  deg;  Alaska-Hawaii  ST,  150  deg. 

The  equation  of  time  can  be  read  from  Table  2  based  upon  the  month 

of  the  year. 

Example:  Calculate  Actual  Solar  time  at  1:30  P.M.,  on  July  21  for 
Manhattan,  Kansas,  longitude  =  96  deg. 

Solution:  Local  standard  time  is  actually  12:30  P.M.  as  Manhattan  is 
in  the  Central  Zone  and  uses  Daylight  Saving  Time.  The 
equation  of  time  from  Table  2   is  -6.2  minutes  and  the 
Central  Zone  is  90  deg  east  of  Greenwich 


AST  =  12:30  +  (-6.2)  -4  (96-90) 
AST  =  12:00  P.M. 
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Table  1. 


Equation  of  Time  for  the  Twenty-first  Day  of  Each  Month 


Equation 

of 

Month         Day 

of  the 

year 

time,  minutes 

January 

21 

-11.2 

February 

52 

-13.9 

March 

80 

-  7.5 

April 

111 

+  1.1 

May 

141 

+  3.3 

June 

173 

-  1.4 

July 

202 

-  6.2 

August 

233 

-  2.4 

September 

265 

+  7.5 

October 

294 

+15.4 

November 

325 

+13.8 

December 

355 

+  1.6 

Source:  ASHRAE  Applications,  1974b. 


Table  2.    Central  Standard  Time  and  Corresponding  Solar  Time 

for  Manhattan,  Kansas 


Month,  date 

CST 

Solar  time 

January  21 

8:00  A.M. 

7:25  A.M. 

February  21 

8:00  A.M. 

7:22  A.M. 

March  21 

8:00  A.M. 

7:29  A.M. 

April  21 

8:00  A.M. 

7:37  A.M. 

May  21 

8:00  A.M. 

DST 

6:39  A.M. 

June  21 

8:00  A.M. 

DST 

6:35  A.M. 

July  21 

8:00  A.M. 

DST 

6:30  A.M. 

August  21 

8:00  A.M. 

DST 

6:34  A.M. 

September  21 

8:00  A.M. 

DST 

6:43  A.M. 

October  21 

8:00  A.M. 

DST 

6:51  A.M. 

November  21 

8:00  A.M. 

7:50  A.M. 

December  21 

8:00  A.M. 

7:38  A.M. 
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APPENDIX  F 
DATA  INPUT  PROCEDURE 

The  input  data  can  be  divided  naturally  into  four  groups.  These 
divisions  as  well  as  the  input  order  and  format  are  shown  in  Table  4.  The 
first  data  set  contains  those  data  which  are  common  to  all  three  of  the 
options  considered  in  this  study.  The  first  three  items,  solar  heat  gain 
factors,  external  temperatures  and  external  relative  humidities,  are  en- 
tered on  an  hourly  basis  for  one  day  of  each  month.  Thus,  one  day  occupies 
2  full  cards.  Four  sets  of  solar  heat  gain  factors  are  required  in  order  to 
accomodate  the  one  horizontal  and  three  vertical  surfaces  considered.  The 
convention  used  here  is  that  the  first  value  read  in  for  any  particular  day 
corresponds  to  1  A.M. 

The  remaining  three  data  sets  are  comprised  of  the  data  specific  to 
any  one  option.  Each  of  these  groups  can  in  turn  be  divided  into  three 
subsets.  For  each  option,  three  types  of  data  must  be  supplied: 

1)  general  plant  characteristics  such  as  wall  surface  areas, 
room  temperatures,  and  personnel  requirements, 

2)  exterior  transmission  parameters  for  the  vertical  surfaces, 

3)  exterior  transmission  parameters  for  the  roofs. 

As  programmed,  each  subroutine  calculates  the  exterior  transmission  load 
for  the  walls  prior  to  that  for  the  roofs.  Consequently,  the  wall  data 
must  be  entered  first,  as  notated  in  Table  4. 

A  system  has  been  devised  whereby  each  room  and  surface  considered 
in  the  current  plant  configuration  is  assigned  an  identifying  number.  The 
primary  purpose  of  this  is  to  define  the  input  order  of  the  general  plant 
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data.  A  hypothetical  plant  is  shown  in  Figure  17  where  the  Roman  numerals 
refer  to  rooms  while  the  Arabic  numerals  serve  to  identify  specific  surfaces. 
The  general  form  of  a  READ  statement  in  the  program  HEATLD  is  where 

READ  (5,  J)  (VALUE(I),  1=1,  NVALS) 

J  refers  to  a  specific  format  statement.  VALUE  represents  some  vector  of 
required  data  while  NVALS  is  a  constant  equal  to  the  number  of  values  re- 
quired. For  example,  to  read  in  the  necessary  surface  areas  for  the  cold 
boning  option,  NVALS  would  become  13  and  the  value  for  surface  1  would  be 
entered  on  the  card  first.  Surface  2,  3,  ...  13  would  follow  sequentially 
according  to  Figure  17. 

However,  for  the  hot -processing  option  I,  NVALS  would  be  9  because  of 
the  differences  in  plant  layout,  i.e.  disregard  room  III  in  Figure  17. 
For  hot -processing  option  II,  one  must  disregard  both  rooms  II  and  III  in 
Figure  17. 
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APPENDIX  G.  COMPUTER  PROGRAM 


COMMON    SHGF(4,  12  ,24)  ,  I  C(  12  ,24)  ,  TE(  12  ,2  4  )  ,  CF.V  (3  )  ,  «ATT  (3  )  ,  US  6  (  3)  , 
1SAF(3),NPE0PL (3 ) , ENTHL P(  12 , 24) ,ENTH1M 3) , PRODI , PKCD2 , 1 CL I , Q ( 10  )  , 
2D(  10)  .CNiUND.A.EPSLfWJcLTAR.HO.CPECFL.NHLJR.MCNTHl  12)  , SPhT  1 , 
3SPHT2,HP(3)  .NCAHC, TOTAL (12) 

DIMENSION    SRI- A I  L3J  ,TfcMP(4)  ,0(13)  ,TEP  FO  13(12,24)  , rfHC6(12,24)  , 
lSHGF0o(4,  12,2  4),CJEFFH3),CUEFF2(6)  .RELHUM  12,24) 

DATA    CCEFFI/.2405,  .^-4a,.01  7o3,3&4.  716o, .2405, .44d, .0137  7,36  2.3  629/ 

CAT  A    CCEFF2/54.63  29,  12  30  I.  69  ,  5.  1692  2  ,2  3.  3924  ,1 12  66.65,  .4o05  7/ 

PEAD(5,i)IDLT, ISTRT 

1  F0RMATU2.I3) 
NH0UR=24/  IDLT*.5 

Q  ****************************************************************** 

C      REAO  IN  EXTtRNAL  TEMPERATURES  ANO  RELATIVE  HUMIDITIES 

Q  **  **************************************************************** 

READ (5, 2) ( (TEMPCEl I, J) , J=l ,24) , 1=1 , 12 ) 
RE AC ( 5,2)  ( (RHOBII, J]  ,J  =  1  ,24)  ,1  =  1 ,i 2) 

Q  **************************************************************  **** 

C  REAO  IN  SOLAR  HEAT  GAIN  FACTORS,  SF  G  f  (  I ,  J  ,K  ) ,  AS  A  FUNCT  ILli    CF 

C  ORIENTATION,  MCNTH,  AND  TIME 

C 

C  1=1,     WESTERN    EXPOSURE 

C  1=2,    NORTHERN    EXPOSURE 

C  1=3,    EASTERN    EXFCSJRE 

C  1=4,    ROOF 

r  ****************************************************************** 

REAO (5, 2) ( ((ShGFLd( I , J ,K ) , <= 1 ,24  ) , J=l, 12)  ,  1=1,4) 

2  FORMAT (12F6.C) 
CO    7    1=1,12 

DO    7    J=l,24 

1NDEX=J*ISTR1-1 

IF( INDEX. GT. 24)     INCEX=  lNoCX-24 

T0(I,J)=TEXPDB(I ,INJcX) 

7  RELHUMI  I,J)=RFOEU,INJEX) 
DO    8    K  =  l,4 

DO    8    1=1,12 

CO    3    J=l,24 

IN0EX=J*ISTRT-1 

IF(  INDEX. GT. 24)     UDEX=  INDEX-24 

8  SHGF(K,I,J)=SFGFC21K,I .INDEX) 
DO    3    1=1,12 

INDEX=0   ° 

DO    3    J=l  ,NHOUR 

T0(  I, J)=TO<  I  ,INDEX*1) 

REL  hUM ( I , J ) =R  EL  HUM ( I , I NDEX  *  1 ) 

3  1NDEX=INDEX*IULT 
CO    4    K=l,4 

CO    4     1=1,12 

INDEX=0 

CO    4    J-=i,NHOUR 

ShGF(K,I  ,  J)=ShGF(K,I , INDEX* 1 ) 

4  INDfcX=INJtX*IOLT 

Q  ********************************************************** ******** 

C      CALCULATE  ENTHALPY  VALUES 

Q        ********************************  *4c  **************************  ****** 

CO  6  1=1,12 
DO  6  J=1,NHCUR 
KK=1 

IF(TO(i, J). LT. 32.)  KK=4 

VAL=COEFF2(KN)-CCEFF2(KK*i J/(TJ(  I,  J  ) +4P9.69 J-C0EFF2  (.<«.*2  )  *ALCG  (T0( 
ll,J)+459.69) 
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COMMON    SHGF(4,12,24)  ,T  C(  12  ,24  )  ,  TE<  12  ,24  )  ,  CFM  I J  ) ,*ATT(3) ,USE(3)  , 
1SAF13),NPECPL  (3),ENThLP(  12,24)  ,ENTHIM  3  )  ,  PRGDl ,  PRCD2  , 1  CLT  ,  Q  {  10  ) , 
20  (  10)  ,C»NO,ND,A,EPSLN,GcLTAR,HO,CPECFi_,NHCUM,MCNTI-H  12J.SPHT1, 
3SPHT2,HP(3)  .NCARC, TuTAL( 12) 

DIMENSION    SRl-Al  13*  ,TtMP(4)  ,U(13)  ,TEf  FD 13  (12,24)  , RHG8 (12 ,24 ) , 
lSHGFDt3(4,  12,24),CJEFFK3),GUEFF2(6)  .KELHUM  12,24) 

DATA    CCEFH/.24J5, .448 , . 01 7o3 ,3e>4. 7 lfco , .24  05 , . 4  4d,  . J 13 7  7 , 862  .3 62  9/ 

DATA    GCEFF2/54.6329, 12 30  I. 69 ,5. 1692  2 ,23. 3924 ,1 1266 .65, .4o05 7/ 

FEA0(5,l)IDLT,ISTRT 

1  F0RMATU2,I3) 
NHOUR=24/IDLT*-.5 

Q  ****************************************************************** 

C      REAO  IN  EXTtRNAL  TEMPERATURES  AND  RELATIVE  HUMIDITIES 

q      A***************************************************************** 

REAJ(5,2)  (  (TEMPCedf  J)  ,  J  =  l  ,2  4  )  ,  1  =  1  ,  12  ) 
REA0( 5,2)  ( (RHDB(I.J)  ,J=1  ,24)  ,1-1,12) 

Q  ********************************************** ********* *********** 

C  READ  IN  SCLAR  HhAT  GAIN  FACTORS,  SFGF(I,J,K),  AS  A  FUNCT IGN  CF 

C  ORIENTATION,  MONTH,  AND  TIME 

0 

C  1=1,  WESTERN  EXPOSURE 

C  1=2,  NORTHERN  EXPOSURE 

C  1=3,  EASTERN  EXPOSURE 

C  1=4,  ROOF 

C  ***********************************  ***************************  **** 

READ  (5,2)  (  (  (ShGFLdU.J  ,K  ) ,  K=  1  ,24  ) ,  J=  1 ,  12  )  ,  1=  1,  4  ) 

2  FORMAT (  12F6.C) 
CO    7     1=1,12 

DO    7    J=l,24 
1NDEX=J+ISTRT-1 

IFUNDEX.GT.24>     INCEX=  INGEX-24 
T0( I, J)=TtMPDB(I ,  INJEX) 

7  RELHUMI I,J)=RF0e(I,INJEX) 
DO    8    K=l,4 

DO    8    1=1,12 

CO    8    J=l,24 

INGEX=J*ISTRT-1 

1FI INDEX. GT. 24)     IMOE X= INDEX-24 

8  SHGF(K,I,J)=SFGFCaiK,I .INDEX) 
00    3    1=1,12 

1ND£X=0   * 

DO    3    J=l ,NHGUR 

T0(  I,  J)=TO(  I  ,IN0EX*1) 

RELhUM  (I, J)=RELhUM( I  .INDEX* 1) 

3  INDEX=INOEX+IDLT 
CO    4    K=l,4 

CO   4    1=1,12 

INDEX=0 

DO    4    J=1,NH0UR 

ShGF( K,I , J)=ShGF(K,i ,  INDEX* 1) 

4  INDtX=lNDtX+IOLT 

C      ****************************************************************** 

C      CALCULATE  ENTHALPY  VALUES 

C     ****************************************************************** 

CO  6  1=1 ,12 
DO  6  J=1,NHCUR 
KK=l 

IF(TO(I,J) .LT.32.)  KK=4 

VAL=C0EFF2(KM-CCEFF2UK*1 )/  (T0(  I ,  J  )  *t!><i  .691-CCEFF2  l<K  +  2  )  *ALGG  (T  0( 
II, J)*459.69) 
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PS=EXP(VAL) 

PV=RELHUM(  I  ,J)*PS/10C. 
XCALC=.6219*PV/(14.toS6-PV) 
KK=1 

IF(T0( I, J) .LT.32.)    KK=5 

ENTHLPU  ,J)-CCEFflUKmO(If  J)+XCALC*UT0U,J)♦459.69>*(CO£Ff:^IClC^■ 
ll)-C0EFFl(^K^■2)  )+CGEFFKKK*3)) 
6    CONTINUE 

KRITE(6,1C)  ,    T^= 

10    FORMAT  (  '  1'  ,T45,  'OUTSIDE    TEMPERATURES    IN    DEGREES    FAhPE.NhE  I  T  '  ,/,  T55  , 

IMMCNJh    1    IS    JANUARY  l»  ,//#T74i  MCMTH't/1 

DO    15    1=1,12 
15    MONTH! I )=I 

WRITE(6,20) (MCMTH( II ,1=1,12) 
20    FORMAT  (  T6,  'HOUR'  ,T18, 11(12, 7X), 12) 

CO    25    J=l,NHCUR 
JJ=I01T*( J-lJf ISTRT 

IFUJ.GT.24)    JJ=JJ-24 
25    WRITE(6,3O)JJ,(T0(I,J),I=l,12) 
30    F0RMATl/T9,I2,6X,12(F4.0,5X) ) 

WRITE(6,35) 
35  FORMAT (//,T54, 'SOLAR  HEAT  GAIN  FACTORS') 

DO  40  1=1,4 

WRITE Ut 45)1 1  (MONTH  J)  , J  =1,12) 

45  FORMAT  I//, T61, ' EXPCS LRE=', 12, /,T74,' NO  NT H',/,T 8, 'HOUR ',T la, 11( 

1I2,7X),I2) 
DO    40    J=l,NhOUR 
JJ=IDLT*( J-1J+ISTRT 
IFUJ.GT.24)    JJ=JJ-24 
40    WRITE  (6,30)  JJ,  (SFGF  (  I ,  K,  J)  ,.<=1 ,  12  ) 
WRITE(6,4t)  IMCMHI  I  )  ,1=1  ,12) 

46  FORMAT 1//T 53, ' ENThALPY    VALUE S ',/, T 74  ,' MONTH'  ,/,T6,'hOUR'  ,r  13, lit  12 
1,7X) ,12) 

CO    47    I=1,NH0UR 
JJ=1DLT*(I-L)+ISTRT 
IF(JJ.CT.24)     JJ=JJ-24 

47  WRITE (6, 48) JJ,(ENThLP( J, 1) ,J=1,12) 

48  FCRMAT(/T9,I2,6X,12(F4.1,3X) ) 

READ  I  5, 49) NCARC, PRCU  I, PR CJ2,SPHT I, SPHT2.CFE0PL 

49  FORMAT  116,  UF6.C) 
FR001=PRCD1*NCA8C 
PRGD2=PROD2*UCARC 

£,  ******  ********  ******  ********************************************** 

C  COLD    BCNING    OPTICN 

Q  ************************************** **************************** 

RE AD  I  5, 5) (SRFA( I) ,1  =  1,13 ) 
READl 5,5) (U( I ) ,1=1,13) 

READ (5, 5) (TEMP( I), 1=1,4) ,T1EX IT ,T2EX  IT , T3EXI T, (ENTHIfU I )  ,1  =  1,3) 
READ (5,5) (CFM( I)  ,1  =  1  ,3), («ATT( I)  ,1=1,3),  (CSEII  )  ,  1=1,3) , 
1(SAF(  I), 1=1,3) 
RE  AD  (5, 5)  «TRMSTSPHT.»T,T3hKD,hrGWTK 
READ( 5,5)(HP( I) ,1=1,3) 
5    FORMAT (13F6.0) 

READ  (5,55)  (NPECPH  I)  ,1=1  ,3) 
55    FORMAT (316) 

CALL    CLOCPT(SRFA,MODfc,TEMP,TlEXIT,T2EXIT,T3E>lT,HTRMS,SPHT«T, 

1TSHRD,HFG«TR,C) 
IF(MCOE.EQ.l)    GO    TO     10C0 

Q  ********************  ******************************************  **** 

C  HOT    BCMNG    OPTICN    1 
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Q  *****************************************  *  ********************  **** 

READ(5,5) (SRPA( [J ,1=1,9) 
READ(5,5)(Ul  I  ), 1=1,9) 

RE  AD (5, 5)  (TEMP (I), 1=1, 2)  ,T1EX IT, T2EX IT , ( cNTH  IN (  1  ), 1  =  1,2) 
REALM  5, 5) (CFM(I)  ,  1  =  1  ,2 ) , ( nAT T( I )  , I =1  ,2 ) , ( USE (  I )  ,  I  =  1 ,2  )  , 
1  (SAF(l),  1=1,2) 
REAG(i>,5)  (HP(  1)  ,1  =  1,2) 
REA0(5,55)(NPC0HL( I) ,1=1,2) 
CALL    GPTKSKf  A, MCJE, U, TEMP , T LEXIT, T2£X  IT  ) 
IFtMCUE.tC  1)    GC    TC    10C0 

Q  ****************************************************************** 

C  HOT    BCMNG    OPTION    2 

Q  ****************************************  ***************  *********** 

REA015,5) (SRFA(  I),  1=1,5) 

READ(5,5)(U(I)  ,1  =  1,5) 

READ! 5, 5) TEMPI  1 ) , T  It  XI T, ENTH IN ( 1 ) 

R£AO(5,5)CFM(  1)  ,WATT(1 ) , USE ( 1 ) , SAP ( 1 ) 

REAO(5,5)HPl 1) 

REAC(5,55)NPECPLll) 

CALL  CPT2(SRFA,MCUE,1,,TEMP,T1EXIT) 

1000  WRITE(6,1001) 

1001  FCRMATl'l*) 
STOP 

END 

SUBROUTINE  CL COPT (SRF A, MODE, TEMP, T IE XI T, T2EXIT , T3c XI T, WTRMS , 

lSPHThT,TSHRt!,hFGV,rR,U) 

COMMON    SHGF( 4, 12,24) ,  IG(  12,2  4) ,  TEt  12,2  4) ,CFM(3) ,^aTT13)  ,LSEl 3)  , 
1SAFI3)  ,NPE0PL(3)  ,  ENThLH  12  ,24  )  ,  ENTh  IN(  3  ) ,  PROD  1 ,  PRO  02, 
2IDLT,B(  lO),DllO),C,NE,NO,A,EPSLNfDELTAR,HC,CPECPL,.',hCJR,MCNTh(  12), 
3SPHTl,SPhT2,hP(3) ,NCAKL, TOTAL  I  12) 

01  MENS  I CN    SPFAI13)  ,TcMP(4)  ,uESJM(12)  ,i,SUM(3,  12  )  ,  STOR(  3,  12)  ,  U(  1  2) 

REAL    INFlT1,INFLT2,INFlT3 

hRlTEl6,l)  (MCNTMl  ),  1=1,  12) 

1  FORMATC  l»  ,T47, 'COLO    3CNING    OPTION' ,// T45 ,' EXTERIOR    TRANSMISSION    ii 
1Y    SURFACE     (BTU/CAY)' ,//T?4, 'MONTH' , / 1 7 , • SCLROE ' , T18, 11(12, 7X), 12) 

Q  ****************************************************************** 

C  CALCULATE  EXTERIOR  TRANSMISSION 

C 

C  hALLS  1,  2,  6,  10,  11 

C  **************************************************************  »*** 

DO  2  1=1,3 
DO  2  J=1,12 

2  QSUM( I,J)=0. 
REAU(5,5).\B,NC,C,A,EPSLN,H0 

5  FCRMAT(2I3,4F10.0) 

REAC(5,10)  (B(  I),I=1,NB) 
READ(5,1C) iDil) ,1=1, NO) 
10  FORMAT teFlO.OJ 
JJ  =  0 
J=l 
C 
C      KJ  IS  A  PARAMETER  CESIGNAT1NG  ThE  ORIENTATION  OF  THE  SURFACE 

C 

KJ=1 

READ(5,10)DELTAR 

DO    25    1=1,5 

LL=I+JJ 

CALL    TFANEX(SRFA(LU  , T £MP ( J ) ,U ( LL ) ,MCOE, QESUMfKJ  ) 

IF(MODE.EU.I)    RETURN 

CO    20    KK=1,12 
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20  CSUMt J,KK)=C£SUf(KK)+3SUMtJ,KK) 
1FU.EC.1)     GO    TO    21 
IM1.EC.2)    GO    TC    22 
IFCX.EQ.3)    GO    TO    23 
IFU.EC.4)    GO    TO    24 

GO    TC    25 

21  KJ  =  2 

GO    TO    25 

22  J=2 
JJ=3 

GO    TO    25 

23  J=3 
JJ=6 

GO    TC    Z5 

24  KJ=3 

25  WRIT  El  6, 2  5) LI,  IGESLMIL  ltL«lt  12  I 

26  FORMAT (T2,« SURFACE       '  ,  12  ,T  15  , 1 1  (  F8  . G  ,  1 X )  ,  F3  .0  )  ***».**. 
c              ******************  ****************************  *************  ******* 

C  fcOOFS    5 1    S i     13 

c  ********************** ***T** *************************** *********** 

READ(5,5)NtJ,N0,C,A,EPSLN,H0 

frEAC(5,10)  (Ed)  ,1=1.  Mi) 

READ(5,U)(0( I) ,1*1, NO) 

REAClb,10)CELlAR 

J=l 

KJ«4 

CO  35  1=5,13,4 

CALL  TKANEXt  SRFAl  I  )  ,  TEMP  (J)  ,U(  I  )  .MCCE*  CtSUN ,  KJ  ) 

IF(MG0E.EG.  1)  RETURN 

CO    30    KK=1,12 
30   6SUMtJ#Ki<>*4sUi(J,KKI*CESUMUM 

J=J*1 

35  ViRITE  (6,26)1,  I  CESJMU)  ,L  =  1,12) 
k«RITfc(6,3£) 

36  FORMAT  1//T5, 'EXTERIOR    TRANSMISSION    TCTALS    (dTU/CAV)  • ) 

DO    37    11=1,3 

37  *RITE(£,33)  11  ,USUM(  II  ,L),L=1,12) 

38  FORMAT  l//T6,'RM    • ,  12  ,T  15  , 1 1 1 Fo.d, IX ) ,F 8.0 ) 

r  ****************************************************************** 

C      CALCULATE  INTERIOR  TRANSMISSION 

Q  ****************************************************************** 

VAL12=U(3J*SRFA(3)*< TEMP (2 J-TEMP ( 1 ) )*24. 
VAL14=U(4)*SKFAl4)*UEMP(4)-TEMP{  1)  )*24. 
VAL24^Ut3)*SRFA(8)*l  TEMP (4 ) -T EMP ( 2 )  )*24. 
VAL23=U(7)*SRFA(7)*( TENP (3)- TEMP ( 2) )*24. 
VAL34=U(12)*SKF*(12)*(TEMP(4)-TEMP<2)) *24  . 
KRITE(6,4l)  VAL12,VAL14,VAL24,VAL23  ,VAL34 
41  FORMAT (//, To,  'INT  TRAN  • , /, T6, 'RM  l->2 '  , T 7 2 ,F 6. C , / , T6, • RM  l->4'  ,T73 
ltFd.0,/,T6,'RM    2->4'  ,T73,F8.0,/,To, 'KM    2->3 • , T 73 ,F a .0, /, To,  •  RM     3-> 

24',T73,F3.C) 
c  ****************************************************************** 

C      CALCULATE  INFILTRATICN  LCSS 

'  **************************************************************  **** 

00   45    1=1,12 

INFLT1=0. 

INFLT2=0. 

INFLT3=0. 

CO    40    J«l,NriOUR 

INFLT1=I.\FLT1*4.5*CFMI  1)  *I  DL  T  *  (  ENTriLPt  I  ,  J  )-ENTHlN(  1)  ) 

INFLT2  =  i:«FLTz*4.5*Ct-fl2)*IDLT*(ENThLF(  I,J)-EMhIM2)) 
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40  INFLT3=INFLT3*4.5*CFM(3)  *I  DLT*  (  ENTHLPI  1,J)-ENTHINI  3)) 
STOft(l,I)=INFLTl 
STOR(2,I)=INFLT2 
ST0R13, 1  )=I.\FLT3 
CSUMI  ltl  )=vjSUM(  1,1  J  +  INFLTl 
CSUM(2,I)=Q3UM(2,I)+I,\FLT2 

45  CSUM(3,1)=GSUM(3,I)+INFLT3 
hRITE(6,44) 

44    FORMAT 1//T4,  '  INFILTRAT  ION'  ) 
00    42    J=l,3 

42  WRITt(6.43)J,( STCRI J, I ), 1*1,12) 

43  FORMAT  i/T8 1  *RH    ',  12.T15,  IK  F8.  L,  IX  ),  F6  .1 ) 

************************************* 4*********44*44**4* ********** 

C      CALCULATE  PERSCNNEL  LOAD 

Q  ********************** ******4*******<*****4*****4*4*44****** 4* ***# 

£PERSl=NPtGPL(  1  )4CPEC?L*8. 
CPERS2=NPECPL (2  )*CPECPL*8. 
gPERS3=NPECPL(3)*»FtCPL*8. 
WRITE(6,4o)CPERSl,CPtKS2,wPERS3 

46  FORMAT  {// ,  T6  ,  •  PtnSCNf^L'  ,/,T3,  «RM    1  •  ,  T  73 ,  F  d  .  C,  /  ,  Td  ,  «RM    2',T73,Fd.G 
2,/,T8,'RM    3',T73,FS.C) 

Q  ****************4*44*********4****4*****4*4***4**4*4 ****** ******** 

C      CALCULATE  PrtCCuCT  LGAJ 

Q  *****  444*  *************************************  4*4**4******  ******** 

PROO=PROJl*SPhTl 
VALl=(Sc.5-TlEXIT)*P»-CC 
VAL2=(TlEXlT-T2EXIT)*PRaO 
VAL3=(T2EXlT-T3£XlT)*P80i:2*SPhT2 
hRITE(6,47) VALl,VAL2,VAL3 

47  FORMAT {/,13,'PRCCJCT    LOAD'  ,/Td,»RM     1 • ,  T73 ,F9 .0, / ,T 3, »RM    2* ,T73  ,F9. 
10,/,T8,'R*    3'  ,T73, F9.0) 

r  ****************************************************************** 

C      CALCULATE  ELECTRICAL  MCTCR  LOAG 

£        *********** **********-****************4*****4*********4*4********* 

ELECT  1  =  HP( 1) *2545. *24. 

ELECT2=hP(2)42545.42*. 

EL£CT3=HP(3)42545. 

WRITE (6, 4d)ELECTi,ELECT2,ELECT3 

48  FORMAT (/,T3.« ELECTRICAL    LOAJ « , / , Td ,  ' KM    1  •  ,  T73, Fd  .0 , /, T o, «RM    2',T73 
l,F8.0,/tT8,»krt    3«,T72,F3.0) 

q  **********************  ********** ******************  ********** ****** 

C      CALCULATE  LIGHTING  LCAC 

Q  **************** **************************** ********** ******  ****** 

ALGHTl=hATTm4LScll)*SAF(  I) *dl.84 
ALGKT2=WATT{2)*USE(2)*SAF(2)*81.e<* 
ALGHT3=WATTi3)*USE(3)*SAF(  3)*ol.84 
hRITE(6,49) ALGhT 1 , Al GHT2 , ALGHT 3 

49  FORMAT!/, T4, 'LIGHTING  LOAD • , / , TO , ' Hf    L» , 1 73.F8.0,/  ,T3, »RM  2',T73,F 
18.0,/,T8,'RM  3,,T73,F8.0) 

r      ***********************************  *  **************  *****  *****  ****** 

C  CALCULATE    SHROUD    LOAD 

£  ****************  ******  ***********Jt  **************************  ****** 

QSHRU=WTRMS*SPHT«T*(  IShKD-TEMPJ  L) ) +HFG wTR4hTPMS 
WRITE16,51 JCSFRQ 
51    FGRMAT(/,T5,'SHRCJU    LOAJ • ,/ , Td , » M    i',T73,F8.0) 
DO    50    I=L,12 

CSUMll  ,I)s(viSUM(  1,1)  ♦VAL12*-VAL14+CSFR0«wP£RS1+VALH-ELECT1+AlGHT1>/ 
1PR0U2 

CSUMI2,I)=IQSUM(2,  D-VAL  12*VAL24*VAL23  +C.PERS2+  VAL2  SELECT 2  + ALGH  T2  )  / 
IPR0J2 
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50    CSUM(3,  I  )=(0SUM<3,  I ) -VAL  23* VAL  34 +  QPERS  3* VAL 3 4-ELECT3 «-ALGHT3  ) /PR  CD2 
WRITE(6,55) ( (CSUMt I , J) , J=l , 12) , 1=  1 ,  3  ),  PRGC  1 ,  PRCJ  1,  PRUJ2 

55  FORMAT(/,T52, '***T0TAL  FCR  COLO  dLM  NG*»*  ■  ,/ 1 1  52  , '  (  2TU/l3  UF  F  IN  IS 
IHEO  PRCCUCTS/T8,  •  RM  1  • ,  T I  5 ,  1  l(  F  6.  1 ,  IX  J ,  F  6.  1 ,  /  T  d  ,  •  RM  2  •  1 115  ,  11  (Fa. 
21tlX),Fa.l,/Ta,,RM3,,T15,ll(Fa.l,lX),Fe.l,//,T/7,  •  fcOUH  '  »/T52,  'l', 
3T73, •2',T104,'3,,/T17, 'PRGJLCT  MASS  (LBS  J •  ,T49  ,F9. 0  ,  17  4  , F9.0  ,T  100, 
4F9.0) 

00    56    1=1,12 
T0TAL( I)  =  0. 
DO    56    J=l  ,3 

56  TOTAH  I)  =  10TAL(  I)*CSLM(J  ,1  ) 
WRITE (6, 57) (TOTAL ( I) ,1=1,12) 

57  FGRMAT(/T2, 'TOTALS    BY    MONTH  •  ,/  T  15  ,  11  (Fa.  1 ,  IX  ) ,  Fa  .1  ) 
WRITE(6,60)NCARC 

60    FORMAT  (/T20 . • NJMBER    LF    CARCASSES    PER    DAY    :       SMI 

RETURN 

END 

SUBROUTINE    OPT  1  (SKFA ,MCDE, U, TEMP , T IE XI T, T 2EXIT ) 

COMMON    SHGF(4,12,24)  ,T0(  12  ,24) ,T£( 12 ,24 ) ,CFM (3 ) , KATT13 ) , JSt ( 3)  , 
1SAF(3),NPE0PL(3),ENTHLP(  12,24)  ,ENTHIN(3)  ,PR001  ,PkuD2 , 1 CLT , d 1 10  ), 
20(10)  ,C,N6,NC,A,LPSL.\,  CElT  Ak  ,  HJ  ,  v.P  EG  PL  ,NhGUR  .MCNTri  I  12)  ,  SPHT1, 
3SPHT2,HP(3) tNCARC, TOTAL! 12) 

0IMENS1CN    SRFA(13),TEMP(4),QSUM(3, 12 ) , QESCMt 12) , STOR ( 3 , 12) ,U ( 1 3) 

RfcAL    INFLTUINFLT2 

WRITE(6,1)(MLNTH( I ), 1=1,12) 

1  FORMAT  ('  I*  ,T47,  •  hOT    BGMNG    OPTION    1  '  ,//T45,  •  tXTERI  Orl    TKANSMISSIGN 
1BY    SURFACE    (BTL/CAY)  • , // T74, • MCNTH •  ,/l 7,  • SCURCE'  ,Tld,ll(I2,7X),I2) 

C  ************  *****&1fZ**$**#ZZZ*t****1l  4*  **#****%  4$  **$*■*******  ******* 

C  CALCULATE    EXTERICR    TRANSMISSION 

C 

C  WALLS    1,    2,    6,     7 

Q  ******  ***$**#*******  ******t.*zrf«*it*$tt>r**>f*4*$Z*Z*Zt£b**  *****?■*£**£ 

DO    2    1=1,3 
DO   2    J=l,12 

2  GSUM(I,J)=0. 
REA0(i>,5)N0,NC,C,A,EPSLN,H0 

5    FCRMAT(2I3,4F10.C) 

REAO( 5,10) (B(  I)  ,1=1, N3) 

READ15.10) (CI I),I=1,NC  ) 
10    FCRMATUCF6.0) 

RCAO( 5, 1010ELTAR 

KJ  =  l 

J»l 

JJ=0 

CO    20     1=1,4 

LL=I+JJ 

CALL    TRANEX(SRFA(LL)  ,  TEMP(  J)  , U(  LL)  ,  f  Ci>£  ,  CESU^  ,  K  J) 

IF(MCDE.Ew.l)     RETURN 

CO    25    KK=1,12 
25    CSUM( J,KK)=QESUM(KK)+QSJM( J,KK) 

IF(I. £(..!>  GC  TC  26 
27 
Id 

GO    TO    20 
26 

27 
28 


IF(I  . 

EC. 2) 

GO 

TO 

IFCI. 

EQ.3) 

GO 

J  j 

GO  TO 

20 

KJ  =  2 

GO  TO 

20 

JJ  =  2 

J=2 

GO  TO 

20 

KJ  =  3 
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20    fcRlTE(6,29)LL,(CESUM(L),L=i,  12) 

29  FORMATU2, 'SURFACE    •  ,1  2,  T15  ,  11  (  F6.  C  ,1X  )  ,F8  .0  ) 

C  ******************** ************************************** ******** 

C  ROOFS    5,    9 

Q  ****************************************************************** 

RE  AD ( 5 , 5 ) NU ,N C, C i A  , E  PS  LH , HO 
REA0(5,10J (D(I) ,1  =  1, Nri  i 
READ{5,10)(0( I), 1=1, NO) 
READCS.lOjLtLTAR 
KJ  =  4 
J=l 

00    30    1=11,12 

CALL  TRANEM  SRFA(I) , TEMP { J ) , U(  I  )  ,MCCE  ,  tEStf  ,  K  J  ) 
DO  35  KK=1,12 
35  CSUMiJ.KK)»CSUN(J»KK)+CESUM(KKl 
J=J  +  l 

30  WRITE (6,29)1, (CESJML) ,L=1,12) 
WRITEl t,33) 

33  FORMAT  1//T4, 'EXTERIOR  TRANSMISSION  TOTALS  IBTL/CAY) »J 
DC  31  1=1,2 

31  WRITE(6,22)I,(GSL.MiI  ,LK)  ,LK=1,12) 

32  FORMAT  1/16, 'RM»,12,T15,11(F8.0,1X),F£.0) 

Q         *******************V** ********* *********************************** 

C      CALCULATE  INFILTRATION  LOSS 

C  a***************************************************************** 

DO    45    1=1,12 
IiMFLTl=0. 
INFLT2=C. 
DO    40    J=l,NhOUR 

INFLT1=INFLT1+4.5*CFM(1)*IENTHLP(I , J J-ENTHlN ( I) ) *I OlT 
40    INFLT2=INFLT2«-4.5*CFM2)*(tNTHLP(I  ,  J  J-cNTFIN  (2  )  )  «I  OLT 
ST0R1 1,I)=INFLT1 
ST0R12,I)=1NFLT2 
QSUM(  1,I)  =  |.SUM{  1,1  J+INFLT1 

45  CSUMI2,  I)  =  CSUM(2,  I  1+HFLT2 
ViRITE(6,44) 

44    FGRMAT(//T4, 'INFILTRATION' ) 
DO   46     1=1,2 

46  WRITEU.4  7H  ,(STCR(l  ,J)  ,J=l,l2) 

47  FORMAT (/T 8, 'RM • , I  2, T  La , 1 1( Fd. J , 1 *)  , F c. 0 ) 

C  CALCULATE    INTERIOR    TRANSMISSION 

£  ****************************************************************** 

VAL12=U(J)*SRFA(3)*(Tw^P(2)-TtMP(l) )*24. 
VAH4=U4)*SRFA(4)*(  TEMP (4)- TEMP ( I ) ) *2  4. 
VAL24=U(d)*SKFAl8)  *(TEMP  14)-TE-1P(2)  )*24. 
WUTC(6,4d)  VAL12, VAL14, VAL24 

48  FORMAT  I//, T6,  • INT  TR A  4  ' , /Tb, 'KM  l->2 ', T 73  ,F3 .0  ,/ T6 ,' SM  1->4',T73,F 
18.0^/T6,*f<^  2->4«  ,r/3,F<J.O) 

Q         ***********************************«**************«***********:**** 

C      CALCULATE  PRCCUCT  LGAJ 

Q         **********.*******************.******  «** **********  **  **********  ****** 

VALl=PRClH*SPHTl*(9d.5-T IE  XI  T) 
VAL2=Pf.OJ2*SPhT2*(TlEXir-T2EXir) 
hRITE(fc,49) VALlf VAL2 

49  FORMAT!/, T2, 'PRCCUCT  LGAD' , / , Td, 'RM  1  •  ,  T  7  2  ,F  9.  C  ,  /  ,  T  6 , '  RJ»  2',T72,F9 
1.0) 

C      ****************************************************************** 

C      CALCULATE  PERSONNEL  LQAJ 

Q  ****************************************************************** 
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CPERS1=NPECPL  ll)*6CG..*8. 
CP ERS2=NPECPL( 2 ) *^ Ft CPL* d. 
WUTEU,5UePERSlfQPE?tS2 
5i    FORMAT!//, T6, 'PERSONNEL', /,T3,  'RM    1  •  ,  r  73  ,  FS.  C,  /  ,  Td  ,  »KH    2',T73,F9.0 

C  ****************************************************************** 

C  CALCULATE    ELECTRICAL    MOTOR    LOAD 

Q  ****************************************************************** 

ELECT1=HP(1)*2545.*24. 
EL£CT2=HP(2)*2545.*24. 
WRITE16,52)ELECT1,ELECT2 

52  FORMAT!/, T2, 'ELECTRICAL    lOAO',/,Tj, 'KM     I'  ,T73,F8.G,/,To,  'RM    2»,T73 
1.F8.0) 

c  ******************************************************  ************ 

C      CALCULATE  LIGhTING  LOAD 

C      ****************************************************************** 

ALGHTi=V»ATT(l)*USE{l)*SAF(  11*31.84 

ALGHT2=hATT(2)*LSE(2 I* SAF (2 ) *dl . 8* 

KRITE(6,53)ALGHT1,ALGHT2 

53  FCRMAT(/,T4, 'LIGHTING  LUAD' , / , Td  ,  '  RM  1 ' ,T  13, Fo .0,/ , Td,  'RM  2',T73,F 
18.0) 

CO    50     1=1,12 

GSUM(  l,I)  =  (USUM(l,  I)  +VAL12*-V/SLl4<-VALl  +  wPERSl  +  ELECr  lfALGFTi  J/PR02 
50    GSUM(2, I)=(gSLM(2, I )-VAL  12+VAL24+  \/AL2«- iPER  S2*t  LECl'2 +AL0HT2  J/PRCJ2 
tiKlTE(6,5u)  ( (QSUM( I, J) ,J  =  1,  12) ,1*1 , 2 ) , PKCL1, PR0U2 

56  FORMAT!/, T51,'TCTAL  FOX  HOT  eCNING  CPTICN  I* ,/ ,T52* » IBTU/L3  CF  FIN 
UShED  PRODUCT  )•,  /T8,  'RM  1 '  ,  T  15  , 1 U F6  .  1  , 1  X  )  ,F8.  I  , /T8  ,  • KM  2',Ti5,lil 
2F8.1.1X)  .Fd.l,//T7  7,'RCOM'  ,/T52  ,  •  1  •  ,T  7ti,  •  2  ' ,  /T  1  7,  •  PRO  JUCT  MASS  (LB 
3S)  ',T49,F5.0,J74»FS.C> 

DO    57    1=1,12 
TOTAL(I)=C. 
DO    57    J»i,2 

57  TOTAL  (I)=TCTAL (  I]  «-QSLM  (J  ,  I  J 
*RITt<6,5d) (TOTAL! I)  ,1=1,12) 

58  F0RMAT(/r2, TOTALS    BY    MJNTH  •  ,/T15,  lllFa.  1 ,  IX)  ,F8.1 ) 
KRITE(6,60).\C/>RC 

60    F0RMATI/T20,  'NUMBER    CF    C/lKCASSES    PER    DAY:       ',14) 
RETURN 
ENO 

SUBROUTINE    0PT2( SRFA ,MCD£, U, TEMP ,T 1  EX  I T ) 

CCMMCN    ShGF!4,12,24) , T C( 12 , 24 ) , TE(  12  ,2  4  )  ,CFM ( 3 ) , VA T T I  3 ) , USE ( 3)  , 
1SAF(3)  ,NPECPL(3),ENThLP(12,24)  ,  cNT  h  I  M  3  ) ,  FHCC  1 ,  PR0C2,  I  OL  T  ,  3  <  10), 

2DI10),C,NB,NJ,A,EPSLN,UELTAR,HO,CPECFL,NnCCR,MClMrH(12)  ,  SPhTl, 
3SPHT2,FP(3) , NL ARC, TOTAL (  12) 

DIMENSION    SRFA  t  13)  ,U(13)  ,TEMP  14  )  ,  JESCM  ( 12  )  ,  >.SUM  (3  ,  12  )  ,5  TOR  (  3  , 1  2) 

REAL     INFLT1   , 

WRITE  (6,1)  !,«CNTH  I  I  1,1=1,12) 

1  FORrtAI('l',r56,'H0T    BCNING    OPTION    2 ■ ,// , T45 , • EXT tR I CR    TRANSMISSION 
I    BY    SURFACE    (  BTU/CAY  ) '  ,//T  74  ,  '  MGiMTh  •  ,/  T7,  •  SOLRoE  » ,  T  Id,  1 11  I  2  ,  7X  ) ,  12 

C  ***************************************** ************************* 

C  CALCULATE    EXTERIOR    TRANSMISSION 

C 

C  WALLS     I,    2,    3 

C  ****************************************************************** 

c  ***************************************************************** 

CO    2    1  =  1,  12 

2  CSUM(1,I)=G. 

REALM  5,5)NB,NC,C,A,EPSLN,H0 
5    FGRMAT(2I3,4F10.0) 
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REA0!5,10)(dt I) ,1=1, Nd) 

REA0(5,10){Cl  1),  1=1, ND) 

REA0!5,10)uELTAR 
10    FORMAT!  1CF8.C) 

KJ  =  1 

DO    ZQ    1=1,3 

CALL    TRANEXISRFAd  ), TEMP  11)  ,Ul  I  )  ,MUDE,  QESLI*  ,  KJ  ) 

IF(MODE.EC.l)     RETJPN 

DO    25    J=l,12 
25    CSUM!  1,J1  =  UESUMU)*QSUM  l,Ji 

KJ=KJ+1 

20  WRITE  U, 21)1,  CttfcSUMC  LI  ,L  =  1  ,12) 

21  FORMAT (T2, 'SURFACE',  12, T  15, 11(F3.G,1X),F6.C) 
READ! 5,51NB,NJ,C,A,EPSLN,H0 
READ(5,10)(d!  I),  I-I.N8 J 

REAJ(5,10l(Dt I ) ,1=1, NJ) 

READ( 5.10JDELTAK 

KJ=4 

Q  ****************************************************************** 

C      ROOF  5 

Q         ****************************************************************** 

CALL  TPANCX!SRFM5),TEMP<1),U15)  ,MCCE,  £ESLM,KJ) 
IF(MGDE.EU.l)  RETURN 
1=5 

WRlTE(6,21)It  UESUML)  ,L=1,12) 
DO  30  1=1,12 
30  CSUM (1,1  )=CSUM!1,I  )+CtSJM( II 
1=1 
WRITE  U,  22)  I,  (CSUM(1,L  ),L=1,  12) 

22  FORMAT!//,  T5  ,  '  EXTER  I  CK.  TRANSMISSION  TOTALS  (BTU/DAY  )  • ,  /T  d,  'RM'  ,  12, 
1T15,  1KF8.C.1X)  ,F3.0) 

r     **********************  *************************************  ******* 

C      CALCULATE  INFILTRATICN  lCSS 

Q  ************************  ****************************************** 

DO  40  1=1,12 

1NFLT1=0. 

CO    35    J=1,NH0UR 
35    INFLTl  =  lNFLTl+H.5*CFMll)*IDLT*(ENTFLFt  I, J  J-ENTHINl  1)) 

STOR(  1,  U=INFLT1 
40    CSUMIl,  I)  =  CSUM(1,  D  +  INFLT1 

1=1 

WRITE (6,23) I, (ST0R(1,J ),J=1, 12) 

23  FORMAT!//  ,T+  ,  «  INF  I  LT PAT  I  ON • , / , T S , • RM • ,  I2,T15,lllFd.0,lX),Fd.C) 

Q  ******************** ***£ *******************************  *********** 

t  CALCULATE     INTERIOR     TRANSM  I  SS  ION 

C  ***********************•***************************************#*** 

VAL14=U(4)*SRFA(4)*!TEI>'P('4)-TE.'1P(  1)  )*24. 
WRITE(6,24  )VAL14 

24  FORMAT!//,  T6,  '  INT    TRAN ' , / , T6 , ' KM    l->4»  ,T 7 2 ,Fe  .0  ) 

£  ****************************************************************** 

C      CALCULATE  PERSCNNSL  LOAJ 

Q  ****************************************************************** 

CPERSl=NPECPL(l )*QPEOPL*d. 
*RITE(6,26)wPEkSl 
26    FORMAT!//, To,  'PtRSGNNEL' ,/,Td,  'RM     I  '  ,T  7  3  ,  F8-.  C  ) 

Q  ********************  ****** ******************************** ******** 

C      CALCULATE  PRODUCT  L040 

£      ****************************************************************** 

VAL1=PR0D2*SPHT2*(Q8.5-T1LXIT) 
WRITE(6,27)  VAL1 
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27    FORMAT!/, T2, • FKCCUCT    LCAC ,/ ,Td , *RM     L' tT 73, F 10.0) 

£  **************************************************************  **** 

C  CALCULATE    ELECTRICAL    .-1UTGR    LOAD 

Q  **************************************************************  **** 

ELECT  1=HPI1)*2 54 5.*24. 
WUTE<6,28)ELECT1 
23    FORMAT!/, T3, 'ELECTRICAL    LOAD*  ,  /  ,  I  3  ,  '  RM    1»,T73,F8.0) 

Q  ************************************  ******************  **  ********** 

C      CALCULATE  LIGHTING  LCAC 

C      **********************  ******************************************** 

ALGHTl  =  «ATTU)*USt!l  )*SAF!  11*81.84 
fcRITfctfc ,31)ALGHT1 
31    FORMAI </,T4, »  LIGHTING    LOAD'  ,  /,  To  ,  '  RM    l',T73,F3.0) 

00    45    1=1,12 
45    CSJrt!  1,I)=!USUM(  l,I)+VAH4  +  UPERSH-VALl*ELtCTH-ALGHTl)/PRCD2 

WRITE  16,29) IGSUM! L, 1), I* It  12 J ,PKC02 
29    FORMAT!/, T51,'  TCTAL    FCR    hCT    BCNINS    CFTICN    2  '  ,  /  ,  T52  ,  '  (3TU.L  3    CF    FIN 
1ISHED    PkGGUCT) •  ,/T8,  »RM    l« ,T lS» 1 1( F 8.1 tlXJ ,F8. 1 , // T77 , ' KCuM'  ,/  179, 
2,1« ,/T17, 'PRODUCT    MASS    (  Li3S)  »  ,  T7  5,  FS  .0  ) 
*RITtt 6,5C)NCAKC 
50    FORMAT  1//T20,  'NUMBER    UF    CARCASSES    PER    DAY  ',14) 

RtTURN 
END 

SUBROUTINE    TRANEX! AREA, TEMP, J, MODE, «ESUM,KJ) 

COMMON    SHGFt-*,I2,24)  ,TG(  12  ,2-t )  ,TE(  12  ,2  4  J  r  Cf  M<3  J  ,  WATTI3)  ,  USEJ  3)  , 
1SAFI3),NPE0PL!  JJ.cNTriLP!  12,24)  ,ENTHlM  3)  ,  PRCDi  ,  PRU02  , 1  OLT  ,  d  ( 10  ), 
20! 10)  ,  C,.\  3,  Nu,  A,  EPSLN,  DELT  AH,  HO  ,  iiP  EC  PL  ,  NHCUR  ,  MONTH  !  12)  ,SPHT1, 
3SPHT2.HPJ3J ,NCARC, TOTAL! 12) 

DIMENSILN    QESUM{12),wE(12,24),TEAVGll^J 
C,  ****************************************************************** 

C  CALCULATE    SCL-AlR    TEMPERATURE 

Q  ****************************************************************** 

DO    5    1=1,12 

CO    5    J=1,WH0UR 

S0LRAu=1.15*SHoF(KJ, I, J) 

TE(1, JJ=A*SCLRAD-(EPSLM=OELTAR)/HO*TG( I ,J) 
5    G£li,J)=0. 

CO    10    1=1,12   . 

QESUM( I)=0. 

JCOUNT-1 

ICGUNT=1 
15    CO    20    J=1,NHCUR 

BSUM=0. 

INDEX=J+1 

CO    25    K=1,NB 

IND£X=INDEX-1 

1F(INDEX.LE.C)INDEX=NHCUK 

fiSUM=aSUM«-L(M*Ttl  I,  INCEXJ 
25  CONTINUE 

DSUM=0. 

INDEX=J 

DO  35  K=2,N0 

INDEX" 1«0£X-1 

IFUNDEX.LE.  J)  INJEX=NHOUR 

DSiJM=DSU.'«*OUJ*QEl  I  ,  INDEX)  /  ARE  A 
35  CONTINUE 

CEII ,J)=A«£A*(BSUM-DSUI^TEMP*C) 
20  CESUM!  I)  =  wESUM{  I )  +CE  { I  ,  J  )  *  iDLT 

IF(  JCCUM.E^.l)  GG  TC  50 

CH£CK=lUfcSUtfi  I  )-CSIQR)  /QSTCR*IOO. 
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IFJCFECK.LT.l.)    GO    TO    10 

IF( ICCUNT.Gfc.20)    GC    70    4C 

50    CST0R=tESUMU  ) 

JCGUNT=2 

IC0UNT=ICCINT*1 

CGSUMt I)=0. 

GO    TC    15 
40    V*kITE(£,45HCCUNT,1 

45    FORMAT  (/// ,5X  , 'NC.MCONV  ERGENCE    AFrtR',13,  IX,  •  ITERAT  IONS  •  ,///,5X,'I= 
1«,I3) 

M00E=1 

RETURN 
10    CONTINUE 

MCDE=0 

RETURN 

END 
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ABSTRACT 

A  prevalent  practice  in  beef  packing  plants  is  to  slaughter 
the  animal  and  chill  the  carcass  for  approximately  24-72  hours.  Fre- 
quently, following  the  cooling  period,  the  excess  fat  and  bone  are 
removed  as  the  carcass  is  fabricated  into  muscles  and  muscle  systems 
before  packaging  and  shipment.  Recent  studies  on  product  quality 
indicate  that  a  procedure  usually  termed  hot  boning  or  hot  processing 
can  be  applied  to  produce  beef  steak,  roast,  and  ground  products. 
Hot  processing  entails  the  removal  of  excess  fat  and  bone  from  the 
beef  carcass  soon  after  slaughter.  The  traditional  cooling  period 
after  slaughter  is  either  shortened  or  eliminated,  thus,  the  carcass 
is  warm  when  cut  and  the  resultant  cuts  are  subsequently  chilled. 

This  study  compares  two  hot-processing  procedures  to  conven- 
tional centralized  beef  processing  to  determine  what  method  is  the 
most  energy  efficient.  A  computer  program  was  written  to  calculate 
the  total  refrigeration  energy  of  each  option.  The  results  showed 
that  a  32  percent  energy  savings  was  achieved  for  carcasses  conditioned 
for  8  hours  before  boning  and  a  42  percent  savings  for  carcasses 
boned  immediately  after  slaughter.  This  comparison  determines  the 
difference  in  the  average  energy  requirements  of  the  various  options 
and  is  not  applicable  to  equipment  sizing. 

The  energy  comparisons  are  based  upon  model  facilities  that 
slaughter  sixty  head  of  cattle  an  hour  for  eight  hours  a  day.  How- 
ever, conclusions  found  for  this  size  hot-processing  facility  should 
not  be  unreservedly  ascribed  to  other  sized  plants. 


